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ABSTRACT 


This  research  is  a  study  in  seepage  analysis  -for  engineered 
structures  built  below  the  groundwater  table;  in  particular  an 
analysis  of  under  drain  systems  -for  such  structures.   General 
seepage  theories  and  design  considerations  a.r<s    discussed  as 
part  of  the  literature  review.   Also  design  concepts  and 
methods  of  analysis  a.r&    covered  which  gives  a  brief  overview  of 
design  considerations  and  techniques  which  are  currently  used 
for  this  type  of  seepage  problem. 

Key  variables  to  the  design  of  an  under drain  system  are    studied 
in  detail  by  using  a  two  dimensional,  steady  state  finite 
element  model..   Such  variables  studied  ar&    drain  spacing,  drain 
width,  anisotropic  hydraulic  conductivity,  thickness  of  the 
aquifer,  head  differential  and  sloping  aquifers. 

Results  are    presented  in  the  form  of  dimesionless  plots  and 
nomographs  for  drain  discharge  and  free  surface  plots  for  the 
different  variables  studied. 
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CHAPTER  ONE 
INTRODUCTION 

Far  civil  engineered  structures  which  interact  with  soil.,  a 
major  consideration  for  design  is  the  presence  of  groundwater. 
Dams,   buildings,  dry  dock  -facilities,  underground  structures, 
pavements   and  earth  retaining  walls  must  be  designed 
considering  the   quantity  of  seepage,  the  position  of  the 
free  surface,  seepage   forces  and  seepage  control.   Good 
drainage  incorporated  in  the    design  of  these  and  other 
structures  can  often  benefit  through    substantial  cost 
reductions  as  well  as  increasing  the  safety  and  usefulness  of 
the  structure. 

Many  embarassing  and  sometimes  disastrous  problems  often  result 
from  man's  endeavors  to  control  seepage.   Saturated  soil 
backfills  have  caused  many  retaining  wall  failures  due  to 
increased  hydrostatic  pressures  following  extensive  rainfall. 
Highways  and  roadways  prematurely  deteriorate  and  fail  due  to 
trapped  water  in  the  pavement,  section.   Basement  slabs,  pond 
liners  and  drydock  facilities  have  been  severely  damaged  by 
hydrostatic  uplift  pressures  and  seepage.   Most  dam  and 
slope  failures  can  be  attributed  to  excessive  seepage  or 
uplift  pressures.   These  examples  and  untold  others  indicate 
that  good  seepage  and  drainage  control  together  with  an    adequate 


understanding  of  groundwater  is  a  key  element  in  the  design  of 
most  civil  engineering  structures. 

Many  civil  engineering  structures  Are    located  in  areas  where 
groundwater  causes  problems  of  uplift  and  seepage.   The  design 
therfore  must  consider  all  problems  associated  with  high 
groundwater  and  what  measures  must  be « taken  to  assure  the 
integrity  of  the  overall  design.   Casagrande < 1965) ,  in  talking 
about  engineering  risks,  concludes  that,  "The  margin  of  safety 
that  we  incorporate  into  our   structures  should  bear  a  direct 
relationship  to  the  magnitude  of  potential  losses,  and  it  must 
also  take  into  account  the  range  of  uncertainty  involved." 
This  important  concept  hopefully  will  be  considered  throughout 
the  presentation  of  facts  and  data  obtained  in  this  paper. 

Prior  to  the  20th  century,  analysis  of  seepage  and  the  control 
of  drainage  for  enginerring  structures  was  done  by  experienced 
engineers  who  based  their  designs  on  trial  and  error  methods 
which  sometimes  were  successful  and  sometimes  not.   Though 
experience  has  proven  to  be  very  neccesary  and  helpful  in  the 
design  of  drainage  systems,  very    useful  tools  employed  today 
are    numerical  methods  such  as  finite  element  and  finite 
difference  techniques  which  require  the  use  of  the  digital 
computer.   Their  problem  solving  capabilities  have  made 
possible  solving  very  complicated  and  complex  groundwater 
probl ems. 


The  purpose  of  this  research  is  to  study  the  design  of 
underdrain  systems  for  large  structures  placed  below  the 
groundwater  table.  In  chapter  two,  basic  seepage  theories  and 
design  requirements,  designs  concepts  and  useful  methods  o-f 
analysis  are    presented  as  a  brief  literature  review.   This 
chapter  covers  critical  aspects  of  designs  that  should  be 
considered  and  contains  a  brief  review  of  methods  used  by 
geotechnical  engineers. 

In  chapter  three,  the  finite  element  program  developed  by  Dr. 
M.  Aral,  of  the  School  of  Civil  Engineering  at  Georgia 
Institute  of  Technology,  is  discussed  describing  how  the 
design  of  building  underdrain  systems  can  be  modeled  using 
this  technique.   This  finite  element  program  made  possible 
studying  effects  of  key  variables  such  as  permeability, 
hydraulic  head,  drain  geometry,  depth  of  groundwater,  and 
anisotropic  soils.   Program  verification  is  performed  using 
another  program,  SEEP84,  modified  by  the  Civil  Engineering 
Department  at  Virginia  Polytechnic  Institute  and  State 
University,  for  the  microcomputer.   Dupuit's  solution  is  also 
used  to  compare  selected  cases  with  the  finite  element 
programs. 


Chapter  -four  presents  data  generated  from  the  -Finite  element 
study  such  as  position  of  the  -free  surface  and  flow  quantity 
for  selected  systems  of  underdrains  suitable  for  use  beneath  a 
foundation.   The  data  is  presented  in  the  form  of  charts  and 
nomographs  to  be  used  as  a  guide  for  underdrain  system  design. 

Chapter  five  summarizes  the  results  of  the  study  and  presents 
appropriate  concl ussions. 


Appendix  A  contains  program  documentation  for  the  finite 
element  program  used  in  this  study  which  details  to  the  user, 
documentation  to  the  use  of  the  program.   Appendix  B  presents 
example  input  and  output  data  generated  by  this  study. 


CHAPTER  2 


LITERATURE  REVIEW 


Seepage  Theories  and  Design  Requirements 


Groundwater 

Soils,  sometimes  described  as  a  porous  media,  is  a  mass  made  up 
of  solid  matter  having  interconnected  voids,  passages,  -Fissures 
or  pore  spaces.   These  pore  spaces  can  be  filled  by  a  liquid 
and  or  a  gaseous  phase.   Figure  2.1  shows  several  typical  soil 
and  rock  indices  which  depicts  how  pore  spaces  can  act  as 
conduits  for  water  or  sub capillary  openings  in  which  water  may 
be  held  by  adhesive  capillary  forces  (Meinzer , 1942) , 

The  soil  mass  forms  a  complex  network  of  irregular  passages 
through  which  the  fluid  may  flow.   A  key  parameter  of 
groundwater  movement  is  the  porosity  of  the  soil.   As  gravity 
acts  on  the  fluid,  it  moves  through  the  pore  spaces  at  varying 
rates  through  each  open  pore  space.   From  a  macro  scale,  this 
flow  has  been  mathematically  described  by  several  authors 
(Bear,  1981;  Terzaghi,  1943;  Cedergren,  1977;  Harr,  1962; 
Bouwer ,  1978; Freeze,  1971). 


Figure  2.1   Diagram  showing  soil  and  rock  interstices.   A,.  Well 
sorted  sand  deposit  having  high  porosity;  B.  Poorly  sorted 
deposit  having  low  porosity;  C.  Particles  that  may  be  porous 
themselves;  D.  Deposit  whose  porosity  may  be  diminished  by 
deposition  of  mineral  matter  in  the  interstices;  E.  Rock  that 
is  porous  by  solution;  F.  Rock  that  is  porous  by  fracturing 
(after  Meinzer,  1942) 


Moisture  content  in  soils  can  vary  with  depth  -From  the  surface, 
as  shown  by  Figure  2..  2.   Most  seepage  problems  considers  the 
zone  of  saturation  which  is  the  zone  where  the  soil  mass  is 
completely  saturated,.   The  terms  groundwater  level,  phreatic 
water  level,,  water  table,  •free  surface  or    groundwater  surface 
s^re    all  synonymous  to  denote  the  level  of  water  where 
atmospheric  pressure  exists. 
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Figure  2.2   Subsurface  groundwater  distribution 
(after  Bear,  1979) 


In  dealing  with  capillary  -fringe  in  fine  grained  soils  the 
height  of  water  can  very  significant.    As  the  pore  diameter- 
decreases,  the  capillary  rise  increases.   Although  the  soil  i 
considered  unsaturated,  the  soils  moisture  content  will  be 
higher  with  depth. 


The  zone  of  capillary  rise  and  -fringe  zone  may  contain  large 
quantities  of  water  for  fine  grained  soils.   Water  removal  by- 
drainage,  due  to  capillary  forces,  may  be  difficult  for  such 
fine  grained  soils.   Figure  2.3  shows  how  the  actual  moisture 
distribution  is  approximated  by  a  step  distribution.   This 
approximation  of  the  capillary  fringe  level  above  the  free 
surface,  hc  ,  is  often  neglected  in  some  seepage  problems. 
Several  emperical  methods,  based  on  grain  size  diameter  and 
porosity  of  soil,  are    available  to  estimate  h c  (Mavis  and 
Tsui,  1939;  Pol  ubari  nova,  1952,  1962;  Silin  Bekchur  i  n ,  1958)., 


Observation 
well 


-p  =  ~rd2 


Relatively  high 
ground  surface 


Moisture  distributions 
in  soil: 


Assumed 

top  of 
capillary 

fringe    \ 


approximate 


Capillary 
frinae 


0.5  1.0 

Degree  of  saturation.  Su. 

Pressure  distnbution 


p  =  yd,  >  0 


Figure  2.3   Approximation  of  phreatic  surface  and  capillary 
■fringe  zone.   (after  Bear,  1979) 


Many  textbooks  on  seepage  and  groundwater  discuss  water  bearing 
and  retention  capibilities  of  the  capillary  -fringe  and 
movement  o-f  moisture  through  this  zone  (Terzaghi,  1943;  Kovacs, 
1981;  Sowers,  1979).   However,  when  this  height  o-f  water  is  very 
small  relative  to  the  total  thickness   of  the  aquifer,  the 
capillary  fringe  is  sometimes  neglected.    The  geotechnical 
engineer  should  be  concerned  with  the  unsaturated— saturated 
flow  condition  when  the  stability  of  the  engineered  structure 
is  to  be  analyzed,  such  as  a  dam. 


Groundwater  analysis  for  seepage  problems  must  consider  the 
quantity  of  water  to  be  drained,  the  soils  resistance  to 
drainage,  and  what  effects  the  drainage  may  have  on  the  soil 
and  the  structure.   Movement  of  groundwater  can  be  categorized 


either  as  seepage;  capillary  rise  or  flow;  percolation;  or 
turbulent  -flow.   For  this  study,  only  movement  o-f  groundwater 
by  seepage  is  considered. 

Design  Requirements 

Structures  built  below  the  watertable  must  be  designed 
considering  the  -Following: 

(1)  structural  ability  to  withstand  the  applied  hydrostatic 
water  pressures, 

(2)  drainage  control  which  short-cuts  the  -flow  or  movement 
o-f  groundwater  -from  getting  to  the  structure, 

(3)  drainage  control  which  is  incorporated  with  the  overall 
design  to  remove  the  groundwater, 

<4)  assure  water  velocities  be  small  enough  to  prevent 
excess  particle  migration  and  erosion. 

Savings  o-f  dollars,  time  and  natural  resources  can  be 
recognised  when  structures  can  be  considered  on  dry,  stable 
material; in  addition  to  thinner  basement  slab  thicknesses  may 
be  used. 


10 


Darcy' s  Law 

In  1856,  Henry  Darcy  investigated  the  -flow  of  water  in  vertical 
homogeneous  sand.   This  work  has  developed  the  most  fundamental 
relationship  which  is  still  used  in  seepage  analysis  today. 
Darcy"  s  Law  states: 

Q  =  k  i  A  (2.  1) 

where  Q  (L-*/T)  is  a  volume  discharge  or  flow  of  a  given 
cross-sectional  areA^     A  (L  )  in  a  given  time.   The  discharge  is 
also  a  function   of  the  hydraulic  gradient,  i  (L/L)  and  a 
resistance  coefficient,  k  (L/'T)  ,  the  permeability  of  the  soil. 
Investigators,  Muskat  (1937),  Taylor  (1949)  and  Leonards 
(1962),  have  concluded  that  Darcy' s  Law  is  valid  where 
velocities  are    low  or  when  flow  is  laminar,  i.e.  Reynolds 
number  less  than  10. 

Hydraulic  Conductivity 

The  soils  resistance  to  flow  or  the  coefficient  of 
porporti onal i ty ,  k,  for  Darcy' s  Law,  is  also  called  the 
hydraulic  conductivity.   Hydraulic  conductivity  is  a 
coefficient  which  expresses  the  ease  or  resistance  the  water 
has  flowing  through  the  soil.   It  depends  on  the  soils  matrix 
and  fluid  properties  of  the  permeant. 
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Many  researchers  have  performed  tests  to  attempt,  to  develops 
empirical  relationships  of  the  soils  properties  to  hydraulic 
conductivity-   Some  of  these  properties  considered  a.re    grain  or 
pore  size  distribution,  shape  o-f  the  grains  or  pores, 
tortuosity,  specific  surface  and  porosity  (Mitchell,  1976s  Kozeny 
1927;  Carman,  1956;  Cedergren,  1977;  Kovacs,  1981).   All  report 
that  careful  field  i nvest i agat i on  and  careful  conductivity 
studies  ^re    necessary  for  each  design.   Figure  2..  4  is  presented 
to  show  some  characteristics  of  various  soil  types  and  most 
important,  methods  for  determination  of  selected  properties 
(Terzaghi  and  Peck,  1948). 

In  the  unsaturated  zone,  the  hydraulic  conductivity,  moisture 
content  and  pressure  head  is  extremely  complex  and  not  well 
understood.   Again,  field  measurement  is  the  best  method  for 
obtaining  the  best  parameters,  such  as  the  use  of  t heist or 
psychrometer  (Kay  and  Low,  1970)  or  gamma  ray  absorption 
techniques  (Ferguson,  1970;   Aral,  1983).   However,  many 
empirical  formulas  have  been  developed  to  relate  the  soil 
parameters  to  a  relative  hydraulic  conductivity  of  the 
un— saturated  zone  (Kovacs,  1981;;  Irmay,  1954;  Chi  Ids  and 
Col  1 i s-George,  1950;  Gardner,  1958;  Brooks  and  Corey,  1966; 
Maslia,  1980;  Reeves  and  Duriguid,  1975). 
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2.2  Design  Concepts 
Seepage  Forces 

The  e-f-fects  o-f  water  on  the  soil  and  the  structure  is 
controlled  by  the  pressures  in  the  water.   In  soil  mechanics, 
the  water  pressure  is  defined  as  neutral  stress.   For  water  not 
moving,  or  steady  state,  the  pressure  can  be  computed  -from  the 
basic  equation  o-f  hydrostatics 


h   =   _P    +   2  <2.2) 


where  h,  energy  head,  is  equal  to  the  position,  2,  plus  the 
pressure  head  or  piezometric  head,  P/^w  •    Typically  in 
seepage  designs,  the  velocity  head  is  very  small  and  may  be 
neglected.   Then  the  pressure  head  plus  elevation  is  equal  to 

head,  h   (Cedergren,  1977). 

For  the  condition  that  water  pressure  is  greater  than  the  total 
load,  say  o-f  the  soil  and  structure,  the  pressures  are    termed 
upli-ft  pressures  and  1  i  qui -faction  can  occur.   Ther-fore  the 
control  o-f  upli-ft  pressures  must  be  a  primary  design 
con si  deration. 
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Critical  Velocity 

According  to  Darcy's  Law,  the  discharge  velocity  of  seepage 
through  soil  is 


v  ,   =   k  x  i 
a 


(2.3) 


Seepage  is  only  occur ing  through  the  pores  and  voids. 
There-fore  the  true  seepage  velocity  should  be  a  function  of  the 
porosity  of  the  soil  mass.   Re-writing  equation  (2.3)  the 
seepage  velocity  can  be  expressed  as 


v 


sn 


k  x  i 


1 


( 2  „  4 ) 


where  i\  e  is  the  effective  porosity. 

In  flow  conditions,  the  saturation  of  soil  is  moving  into  or 

away  from  the  unsaturated  soil.   The  capillary  head,  hc  ,  is 

working  with  gravity  to  increase  the  hydraulic  gradient.  For- 
th is  condition  the  hydraulic  gradient  can  be  expressed 


h   +   h 


L 


(2.5) 


Figure  2.5  depicts  this  condition  for  flow  through  a  sand. 
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Reservoir 


Reservoir 


Screen 


i  ss  (h  +  /ic)/L 


■o 

CO 


I  +>^  Reference  line 


I- 


(6) 


Figure  2.5   Force  of  saturation. 

Boiling  occurs  when  seepage  pressures  i n  an  upward  direction 
exceeds  the  downward  -force  of  the  soil.   This  must  be  analysed 
for  all  seepage  designs.   Critical  hydraulic  gradient,  i   ,  can 
be  expressed 


\  T  —     y 


w 


t 


w 


If 


* 


w 


.6) 


where 


s 


W 


total  unit  weight  of  soil 
unit  weight  of  water 


b   =   buoyant  unit  weight  of  soil (NAVFAC, DM-7. 1 ) 
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Fi 1 ters 

A  very  important  design  consideration  is  the  possibility  of 
seepage  forces  moving  erodible  soils  and  rock,  termed  piping. 
Betram  (1940),  under  Terzaghi  and  Casagrande,  established  the 
following  criteria  for  filter  design,,  based  on  grain  size. 


D    (of  filter)  D    (of  filter) 

- <  4  to  5  <    — 

D  „_  (of  soil)  D  ._  (of  soil)     (2.7) 


85 
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This  relationship  is  represented  by  Figure  2.6  to  illustrate 
the  prevention  of  piping  by  filters  (Cedergren,  1977). 


LEGEND 
=  m-place  soil 


•       =  Dttr,  soil  particle, 
entrapped  in  filter 


=  soil  which  has 
migrated  into  filter 
and  is  held  by  Due, 
size  soil  particles 


Nominal  boundary 

before  stabilization 

under  seepage 


r>85  (soil) 


Dy,  (filter) 


Figure    2.6 


Illustration    of    prevention    of    piping    by    filters 
(after    Rankilor,    1981) 
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The  ratio  o-f  D15  o-f  a  -filter  to  the   D  g^   o-f  a  soil  is 
considered  the  piping  ratio.   Many  authors  and  organizations 
have  established  other  criteria  -for  -filters,  similar  to  that  o-f 
equation  (2.7)    (U.S.  Army  Corps  0+  Engineers,  1941,  1955;  U.S 
Bureau  o-f  Reclamation,  1955;  Sherard , etal ,  1963;  Cedergren, 
1977;  Tham kachalam  and  Sakthi vadi vel ,  1974). 

Filters  must  be  designed  to  meet  two  major  requirements;   (1) 
they  must  be  sa-fe  with  respect  to  erosion  and  piping,  and  (2) 
must  have  sufficient  discharge  capi bi 1 1 i ti es  to  remove  the 
seepage  quickly  and  e-f  -f  i  ci  ent  1  y.   This  must  be  done  without 
inducing  high  seepage  -forces  or  hydrostatic  pressures.   The 
right  side  o-f  equation  (2.7)  states  that  the  D  jr  size  o-f  a 
-filter  should  be  at  least  4  to  5  times  the  D-ir  o-f  the  soil. 
This  insures  that  the  hydraulic  conductivity  o-f  the  -filter  and 
drains  is  su-f-ficient  to  prevent  the  build-up  o-f  large  seepage 
-forces  and  hydrostatic  pressures. 

There  is  a  wide  range  o-f  engineering  properties  o-f  materials  to 
be  considered  -for  -filter  drainage  systems  designs.   Properties 
which  may  influence  hydraulic  conductivity  are    grain  size, 
density,  mi  neral  ogi  cal  composition,  the  nature  o-f  the  permeant  or 
the  degree  o-f  saturation.   Much  work  has  gone  toward  methods 
-for  determination  and  -factors  which  influence  hydraulic 
conductivity  (Milligan,  1975;  Lambeand  Whitman,  1969; 
Yemington,  1963;  Bouwer ,  1978).   There  is  still  much  research 
necessary  to  understand  better  this  property  o-f  soil. 
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Trench  Drains 

For  the  application  of  drainage  -for  structures  built  below  the 
water  table,  the  use  of  trench  drains,  which  acts  as  collectors 
of  the  seepage  or  may  contain  collector  pipes,  is  a  very 
efficient  method.   Figure  2.7  shows  the  typical  cross-section 
which  may  be  considered  for  the  application  of  this  type  of 
drainage  system.   Trench  drains  will  be  discussed  further  in 
the  report  and  their  utilization  for  this  application  will  be 
studied  in  more  detail. 


FREE  SURFACE 


TRENCH  DRAINS 


Figure  2.7   Typical  application  of  a  trench  drain  system 
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Geosyntheti  C5 

Today,  new  synthetic  membranes  are    used  to  serve  the  purpose  of 
either  -filtration  or  seperation.   This  two  -factors  Are    very 
important  parameters  -for  most  seepage  problems.   Figure  2.8 
shows  comparisons  o-f  some  of  the  -features  o-f  membrane-wrapped 
drains  to  conventional  granular  pipe  drains.   Also  Figure  2.9 
shows  the  bridging  effect  caused  by  good  seperation  of  the 
larger  materials  and  smaller  particles. 


firannlar  rlrain 


Membrane-wrapped  drain 


*+    *t 


fejfcjLi 


Pipe  drains 
need  frequent 
rodding  points 
Expensive. 


Pipe  must  be 
levelled  by 
engineer,  to 
be  an  exact 

gradient.  This  is  expen- 
sive. Trench  must  be  left 
open  for  this,  and  can  get 
silted  up  during  storms. 

Well-point 
pumping 
may  be 
necessary  to 
lay  the  pipe. 


Sand  and  stone  must  be 
supplied  and  laid  to  a  strict 
specification.  Very  expensive. 

Pipe  is  expensive  to  buy,  and 
must  be  laid  by  hand.  Trench 
must  therefore  be  wide,  and 
timbered  if  deep. 


Wide  trench 
means  more 
excavation  and 
more  backfill  than  is 
technically  necessary 


The  deep  transporter  means 
that  settlement  can  take  place 
without  the  drain  blocking  — 
unlike  a  pipe 

Can  be 
taken  deeper 
because  sides 
do  not  need 
support 


t  *»• 


Cost  of  shoring 
up  open  trench 
is  eliminated 

Granular  fill  only 
supports  trench, 
so  can  be  cheap 
specification 

This  drain  is  very 
much  quicker  to 
construct  &  can 
be  laid  in  flooded 
ground. 

Narrower 

width  because  no  pipe  to  be 
laid,  so  men  do  not  have  to 
go  into  the  trench. 


Figure    2.S       Conventional     granular    pipe    trench    drain    and 
membrane-wrapped    trench    drain.       (after    Rankilor,    1981) 
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Magnified  cross  section  of  permeable  membrane 

filter  drain 


Original 

soil 

structure 


Filter 
zone 
in  soil 
(sometimes 
called 
Filter  Cake) 


Bridging 
network 
of  larger 
particles 


Drain 


Membrane 


Figure    2.9       Bridging    e-f-fects    caused 
material-       (after    Rankilor,    1981) 


by    use    o-f    membrane    filter 
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A  good  representation  o-f  how  geosynthetics  combine  action  of 
seperation  and  -filtration,  is  shown  by  Figure  2.10.   As  water 
has  more  area  of  -flow  with  the  membrane,  the  seepage  can  be 
considered  to  be  at  a  -faster  rate  or  drain  quicker  (Rankilor, 
1981) 


Without  membrane 


Water  tries  to 
escape  through 
small  area 


SECTION 


Clay 
soil 


/     Without  membrane     [ 


With  membrane 


Small  surface  area 
for  moisture  loss 


Large  stone/clay 
contact  area 


Water 

escapes 

rapidly 


SECTION" 


Clay 
soil 


o  o PLAN  o  o 

Clay 


tx 


0 


>1 

/w 


Large  surface  area 
for  moisture  loss 


With  membrane 


/^^ 


Small 
stone/clay 
contact  area 


Figure    2.10       Seperation    effects    by    using    a    membrane. 

(after   Rankilor,    1981) 
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2.3  Methods  o-f  Analysis 

In  analysing  seepage  problems;  geometry,  boundary  conditions, 
soil  parameters,  -Flow  quantities  and  location  of  the  free 
surface  must  be  quantified  by  some  means.   Though  this  involves 
many  considerations,  the  engineer  must  gain  reasonable  safe 
results  for  the  design., 

Through  the  history  of  seepage  analysis  many  techniques  and 
methods  have  been  developed,  such  as  graphical,  analytical  and 
numerical  solutions..  Still  used  today  is  the  graphical  method 
of  flow  net  construction.  However  with  the  complex  geometries 
envoi ved  and  via  the  speed  of  computers,  numerical  models  have 
been  developed  as  very  useful  tools  for  the  engineer. 

Flow  Nets 

Many  researchers  have  developed  methods  for  analysing  flow 
through  porous  media.   The  oldest,  developed  by  Casagrande  and 
Terzaghi  (1937),  is  the  use  of  flow  nets,  a  graphical  approach 
to  seepage  analysis. 
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The  continuity  equation  in  the  form 


Su     3v     9w  (2.7) 

t  + 


9x     9y     dz 

shows  -for  a  Cartesian  coordinate  system,  the  quantity  o-f  water 
entering  an  element  must  equal  that  leaving.   The  terms  u,  v 
and  w  are  discharge  velocity  components  o-f  directions  x ,  y  and 
z»    respectfully. 

According  to  Darcy's  Law,  expressed  by  equation  (2.3),  the 
component  o-f  discharge  velocity  can  be  expressed 


3h  8h  9,h 

u  -  -k   ,     v  =  -k  ,    w  =  -k  (2.8) 

9  x  9  y  9  - 


Substituting  into  equation  (2.7) 


-k(9h/9x)      -k(9h/9v)       ~kOh/9z)  (2.9) 

a  +  a +    3 =  o 

9  x  9  v  d  z 
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If  k  is  constant  for  the  element  and  then  expressed  in  a 
typical  two-dimensional  form,  equation  (2.9)  becomes  the  common 
Laplace  equation. 


2 
9zh 


3x 


+ 


82  h 


9Y 


(2.10) 


=  0 


The  general  conditon  of  flow  analysis  using  flow  nets  is 
represented  by  Figure  2.11.  Rules  and  good  examples  are 
presented  in  many  textbooks  (Cedergren,  1977;  Harr,  1967). 


Equipotential 
lines 


Free  water  surface 
(uppermost  line  of  seepage) 


Figure  2.11   General  condition  for  flow  lines  and 
equipotential  lines  in  flow  net  construction. 
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Pup nit  Assumptions 

For  uncon-fined  seepage  problems   the  engineer  must  calculate 
the  location  of  the  -free  surface.   Many  problems  of  flow  have 
been  approximated  based  on  the  Pupuit  theory  of  uncon-fined 
flow.   Pupuit  (1893)  assumed;   (1)  for  small  inclinations  of 
the  line  of  seepage,  the  streamlines  can  be  taken  as  horizontal 
and  equal  potential  lines  vertical,  and  (2)  the  hydraulic 
gradient  was  equal  to  the  slope  of  the  free  surface  and 
invarient  with  depth  (Harr,  1962). 

Harr  (1962)  gives  many  example  problems  and  has  derived  the 
condition  of  free  surface  elevations  as  shown  by  Figure  2.12. 


//// 


/////. 


////////////////I// //// 


Figure  2.12   Simple  free  surface  notation. 
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In    the    case    of    Figure    2.12    the    Dupuit    -formula    is    expressed    as 

2  2 

h  -       h 

q    =    k         — (2.11) 

2L 

Many  other  reports  and  textbooks  have  presented  charts  and 
nomagraphs  based  on  these  assumptions.   Later  in  this  study  a 
nomograph  for  a  simple  trench  drain  with  blanket  drain 
solution  developed  by  the  Federal  Highway  Administration  (1980) 
will  be  presented. 

Numerical  Analysis 

Flow  through  porous  media  has  been  approximated  by  many  authors 
using  ordinary  non-linear  differential  equations  and  partial 
differential  equations.   Coupled  with  the  speed  o-f  the  computer, 
numerical  solutions  to  many  seepage  problems  have  been 
developed.   It  should  be  pointed  out  that  complex  mathematical 
problems  solved  by  numerical  methods,  Are    only  as  Accurate    as 
the  math  model  can  represent  the  conceptions  and  phenomina  of 
seepage  flow.   In  most  cases,  simplification  is  the  rule,  and 
numerical  models  can  only  be  treated  as  approximations  to  the 
actual  occur ances.   Nevertheless,  the  number  crunching 
capibi 1 i ties  and  time  savings  involved  gives  the  engineer  much 
needed  time  for  the  true  analysis  of  the  seepage  problem. 

The  analytical  and  graphical  methods  mentioned  previously  Are 
sometimes  dif-ficult  and  limited  to  flow  systems  which  are 
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relativly  simple™   In  actual  -field  conditions,  the  soils  may  be 
heterogeneous  and  anisotropic  and  boundary  conditions  very 
complex.   These  conditions  make  closed  -form  methods  very 
difficult  and  many  simplifications  would  have  to  be  made. 
However,  with  the  numerical  models,  these  complexities  can  be 
easily  handled  and  many  more  check  cases  can  be  made  in  much 
less  time. 

Finite  Difference  Methods 

Many  numerical  computer  models  have  been  developed  using  finite 
difference  methodology.   Since  the  classic  work  of  Peaceman  and 
Rachford  (1955),  the  oil  industry  has  sponsored  many  studies  on 
model  development.   A  good  historical  background  of  the  finite 
difference  methods  developed  in  seepage  analysis  is  giving  by 
Huyakorn  (1983).   Worth  mentioning  are    the  works  of  Cooley 
(1971)  and  Freeze  (1971)  in  the  development  of  other  finite 
difference  models.   These  models  solve  for  the  free  surface  and 
treat  the  saturated-unsaturated  flow  as  an  immobile  air  phase. 

Finite  Element 

The  finite  element  method  for  solving  subsurface  flow  has  been 
widely  used  and  researched  be  soils  engineers  and  groundwater 
hydrologists.   The  program  used  in  this  research  and  discussed 
in  more  detail  in  Chapter  3  is  a  finite  element  model. 
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Most  -Finite  element  models  are    built  on  the  Galerkin  method. 
Some  earlier  finite  element  models  -for  steady  state  and 
saturated  flow  were  developed  by  Desai  and  Abel  (1972),  Conner 
and  Brebbia  (1977)  and  Zienkiewicz  (1977).   Pinder  and  Frind 
(1972)  used  isoparametric  quadrilateral  elements  with  the 
Gal er kin's  technique.   Neuman  (1973,  1974)  used  a  -four  node 
quadrilateral  element.   A  good  reference  for  other  works  in  the 
finite  element  models  that  have  been  developed  is  reported  in 
Huyakorn  and  Pinder  (1983)  and  Heijde  (1985). 

Some  significant  works  were  those  of  Taylor  and  Brown  (1967) 
where  a  finite  element  model  was  developed  for  seepage  through 
a  dam  solving  for  the  free  surface.   The  flow  was  only 
considerd  for  the  saturated  zone  in  this  model.   Freeze  (1971) 
demonstrated  that  the  conditions  of  the  unsaturated  zone  may 
strongly  influence  the  position  of  the  free  surface  and  results 
obtained  not  considering  the  capillary  action  for  this  case 
could  be  in  error.   Fredlund  and  Morgenstern  ( 1976)  developed 
improved  constitutive  equations  for  an  unsaturated  soil.   Lam 
and  Fredlund,  etal  (1984)  presented  additional  theories  and 
numerical  development  for  further  modeling  techniques  to 
consider  the  unsaturated  zone  in  seepage  analysis. 
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CHAPTER  3 

FINITE  ELEMENT  STUDY 

3.1  Finite  Element  Program  Theory 

Introducti  on 

Finite  element  analysis  o-f  seepage  problems  can  be  described  as 
•follows  (Huyakon  and  Finder,  1983): 

(1)  A  discrete  number  o-f  interconnected  nodes  to  de-fine  a 
series  o-f  sub-regions  or  elements.   These  can  take  the  -form  o-f 
triangles,  quadrilateral  or  rectangular  shapes  which  are  so 
spaced  to  de-fine  the  geometric  boundaries  and  soil  types. 

(2)  A  matrix  expression  is  developed  to  relate  -flow 
variables  to  each  node  o-f  the  elements.   The  mathematical 
expression  may  be  o-f  the  -form  o-f  either  a  variational  or 
weighted  -function.   A  weighted  residual  -function  is  used  in  the 
Galerkin  method  (Conner  and  Brebbia,  1977). 

(3)  A  set  o-f  algebraic  equations  that  describe  the  entire 
global  system  is  applied  to  assemble  the  element  matrices. 
These  equations  are    derived  to  consider  compatible  conditions 
at  each  node  shared  be  each  adjoining  element. 

(4)  Boundary  conditions  o-f  -flow  are    incorporated  to  reflect 
appropriate  limits  to  the  global  matrix  equations. 

(5)  The  resulting  set  o-f  simultaneous  algebraic  equations 

are    solved  by  some  alogrithm  method.   Many  algorithyms  have 
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been  used  to  ef  f  ici  ent  1  y  take  advantage  of  the  banded  and 
symmetric  -features  of  the  coef  f  i  ci  ent  matrix  (Bathe  and  Wilson, 
1976) . 

The  -finite  element  method  has  great  potential  in  seepage 
analysis  problems,  in  that  the  complex  geometries  o-f  flow 
regions  can  be  easily  constructed  via  nodes  and  elements.   This 
chapter  will  present  the  general  theories  and  -formulation  of 
equations  for  the  model  used  in  this  study. 

Program  Theory 

A  finite  element  program  developed  by  Dr.  Mustus-fa  Aral  of  the 
School  o-f  Civil  Engineering,  Georgia  Institute  o-f  Technology, 
was  used  to  study  the  -flow  and  variables  involved  for 
subsurface  drains  as  earlier  described.   This  program  was  used 
to  see  effects  of  different  parameters  in  the  analysis  performed. 
The  program  was  operated  on  the  Georgia  Institute  of  Technology 
main-frame  CDC  CYBER  computer.   A  complete  user's  guide  for  the 
program  is  included  as  Appendix  A-l.   Also  example  input  and  output 
used  in  this  study  is  included  by  Appendix  B-l.   The  formulation 
for  the  program  presented  here  was  summarized  from  a  previous  study 
(Aral,  Sturm  and  Fulford,  1981). 

This  model  only  considers  the  flow  condition  of  the  free 
surface  being  the  point  where  the  water  level  is  at  atmospheric 
pressure.   In  other  words,  this  model  does  not  treat  the 

31 


unsaturated  zone  of  the  soil.   For  this  study  the  -flow  continum 
is  modeled  to  determine  the  quantity  o-f  flow  at  each  drain  and 
to  determine  the  -free  surface  location  between  each  o-f  the 
drains.   This  model  also  computed  the  velocities  at  each 
element  -for  a  good  check. 

Equations  o-f  Flow 

Darcy's  Law  is  used  as  the  basic  relationship  governing 
groundwater  flow.   Flow  is  assumed  to  be  laminar  and  the 
inertial  forces  are  negligible  compared  to  viscous  forces. 
Rewriting  Darcy's  Law  in  the  cartesian  coordinate  system: 


v   =  -  K  V      h  (3.  1) 


where  V  is  seepage  velocity,  h  is  piezometric  head  and  K  the 
hydraulic  conductivity.  In  the  two  dimensional  case,  equation 
(3.1)  can  be  written  as 


v   =  -  Ki±  ,    i=  1,2     (3.2) 

3x 


where  K.,   is  a  conductivity  tensor. 
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The  continuity  equation  -For  the  flow  in  an  element  of  soil  can 
be  expressed 

(Pv.)    =  ~  (H  P  )  ,   1=1,2   (3.3) 

9x  at 


where  T]       is  porosity  of  the  soil  and   p  is  the  density  of  the 
•fluid.   Assuming  the  fluid  and  soil  particles  are 
incompressible,  equation  (3.2)  and  (3.3)  can  be  used  tor  the 
■following  expression  in  the  two-dimensional  -form: 


3        9h 

*Li  _     =   ° 

3x  i      9  x  i 


i  =  1,2 


(3.5) 


Figure  3.1  can  be  refered  to  for  notation  of  the  key  variables 
used  is  this  study. 
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Boundary  Conditions 

For  -free  surface  problems,  which  is  the  case  for  this  study, 
the  exact  location  of  the  phreatic  surface  is  not  known  at  the 
beginning  of  the  solution.   Location  of  the  free  surface  is  an 
iterative  process  using  two  independent  boundary  conditions 
available  at  the  free  surface-   The  line  of  seepage  at  the  free 
surface  is  the  upper  streamline  of  the  flow  domain.   At  every 
point  along  this  streamline  the  pore  pressure  is  at  atmospheric 
pressure  and  therefore  the  piezometric  head  equals  the  elevation 
head.   Also  the  flow  across  this  boundary  is  zero  because  of  it 
being  a  streamline. 

Considering  the  seepage  domain  as  depicted  by  Figure  3.1,  the 
steady  state  boundary  conditions  can  be  described  as  fallows: 

(a)  Impervious  Boundary  (AB) :   The  program  considers  this 
boundary  as  a  streamline  where  velocity  to  this  boundary  equals 
zero. 

(b)  Constant  Head  Boundaries  (AC)  and  (EF) :   Constant 
hydrostatic  pressure  is  exerted  on  these  boundaries  where  total 
piezometric  head  is  constant,  i.e.  equal  to  the  water  level  in 
a  piezometric  tube. 

<c)  Phreatic  Surface  (CE) ,  <FG)  ,  (HI),  (JK)  and  (LM)  :   This 
boundary  is  a  combination  of  both  cases  (a)  and  (b)  above.   By- 
iteration,  described  later,  both  conditions  must  be  satisfied 
independentl y . 
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(d)  Seepage  Faces  (EF) ,  (GH) ,  ( IJ ) ,  (KL)  and  (MN) :  Seepage 
-forces  occur  at  these  boundary  -faces  where  water  exits  through 
the  soil  continum,  i.e.  trench  drains.   Ther-fore  streamlines  oi 
flow  must  be  normal  to  these  boundaries. 
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3.2  Finite  Element  Formulation 
General  Equations 

The  -first  step  in  solving  the  seepage  problem  by  the  finite 
element  method  is  to  discretize  the  -flow  region  into  a  finite 
number  of  subregions  called  elements.   Using  finite  elements, 
equilibrium  equations  must  be  made  to  handle  the  geometries 
involved.   This  results  in  a  set  of  simultaneous  algebraic 
equations  rather  than  differential  equations.   Figure  3.2 
depicts  the  quadrilateral  elements  used  for  this  program  and 
shows  how  four  si«  nodal  triangles  ^re    formed  internal  to  the 
program.   A  more  detailed  description  is  of  finite  element 
equations  and  formulations  can  be  found  in  Zienkiewicz  (1971), 
Desai  (1972)  and  Aral  (1974,  1976). 


ELEMENT  CORNER 
NODES 


MID  NODES 


Figure  3.2   Typical  element  and  nodes 
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A  finite  element  approximation  to  equation  (3.5)  can  be  written 
through  a  variational  approach.   Assuming  Dirichlet  and  Neuman 
boundary  conditions,  and  considering  the  two  dimensional  case 
of  steady  -flow,  the  -following  equations  can  be  expressed: 


h   =    h 


(3.7) 


and    V 


n 


K  1 1 


8h 


3x 


anl  +  K22 


dh 


9x 


a  J  =  v 
nz/     n 


on  S2  (3.7) 


where  h  is  the  specified  Direchlet  boundary  condition  in  terms 
o-f  piezometric  head;   V    is  the  velocity  in  the  direction 
normal  to  the  boundary;  «nl  and  d-n2     Are    the  direction  cosines  of 
the  normal  with  respect  to  Xj  and  ^   ;  and  Sj  and  S2  indicate 
respective  boundaries  of  the  domain.   From  equations  (3.5)  and 
(3.7),  the  following  variational  statement  can  be  written: 


K 


11 


8h 


9  * 


/ 


8x, 


n 


■22 


9  v., 


Shdxj dx2 


-  V  {   6hdS 


n 


(3.8) 
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Integrating  by  parts  one  obtains 


9h   9  <$h        9  h   9  £h 


+  Kyy f    dx  ,  dx  ~ 


9x  i   9  >; ,       8  x  „   9  ;• 


Vn  6  hdS    =    0  (3.9) 


S2 


Since  -from  earlier  discussed  boundary  conditions  as  -For  this 
problem  studied,   Vn  -  0.   Thus,  equation  (3.9)  can  be  further 
simplified  by  eliminating  the  boundary  intregrals. 

Assuming  equation  (3.9)  is  written  -For  a  single  element,  matrix 
equations  can  be  generated  using  an  interpolation  -function  for 
approximating  h  over  an  element.   The  interpolation  -function 
used  in  this  model  is  given  as 

h    =   Nk   h  k      k    =    1  ,  2 , 3 ,  .  .  .  6   (3.10) 

where  h^'s  are  the  unknown  nodal  values  of  the  dependent 
variable  and  N k  is  the  interpolating  polynomial  used  to 
approximate  h  in  an  element. 
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For  a  typical  six-nodal  triangle  element  these  polynomials  are 
given  as 


N  .    =   £,  (2£,  -1) 

1  1     x 


i  -  1,2,3 


(3.  11) 


and 


Ni+3   "   4Vj 


J  =  2,3,1 


(3. 12) 


where     £, ,    £„    ,     and    £„     are    area    cordinates    defined    by    the 
•following    relations 


?1      "    -  ?   2      "     -  ?3=    — 

A  A  A 


where  Ai  ,,  A2  and  A  3  are  the  areas  of  the  three  subtriangles 
subtended  by  the  point  <P)  and  the  corners,  the  index  on  (A) 
designating  the  opposite  corner  numbers,  and  A  is  the  area    of 
the  triangle  (123)  shown  in  Figure  3.3. 


x«  ♦ 


x\ 

*(",», 

1) 

side  1,  £1=  0 

side  2,  E,2=  (, 

\\  A\ 

A2Ni^ 

e 

SP\ 

l(0,0,l)£-— 

A3 

^  *  f«,».°) 

\  side 

2, 

Ca=     0 

Figure  3.3   Area  coordinates, 
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To  simplify  the  resultant  equations,  the  origin  of  the  local 
axis  which  de-fines  the  local  coordinates  of  the  element,  nodes 
are  placed  at  the  centroid  of  the  element  and  the  principle 
axis  are  inclined  angle  ©  in  the  direction  of  local 
ani  sot ropy. 

Substituting  equation  (3.10)  into  equation  (3.9)  the  final  form 
of  the  approximating  equations  in  an  element  can  be  expressed. 
The  stationary  condition  of  equation  (3.9)  is  obtained  by 
equating  the  first  partial  derivatives  of  the  function,  with 
respect  to  the  undetermined  parameters,  h.»  ,  to  zero: 


9Ni   8Nj        cN±       3Nj 


I e  ( h  )  ,  h.     =  \ \    SK,  K22 -|  h±  d x  x  d ^   , 


3*1   9«1        3:=2   9:'2j 


i,j  =  1,2,3,  ..  .6  (3.  14) 


This  results  in  a  six  by  six  element  stiffness  matrix,  S  ,  and 

e 

a  six  by  one  element  load  vector  F. 
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The  six  by  six  element  matrix  equations  of  equation  (3.14)  is 
wr  i  tten 

CSe]    Che>   ■   CFe3-  (3.15) 

where 

9  N  3N.  9  N.      9  N.J 

rfl     -     \\      |K„     — 1      —  J+        K22— 1     — J  dMjdMj, 

9  X  .   9  X  .  9  X  2   9  X  _  J 

i,  j  =   1,2,3, .. .6  (3. 16) 

The  same  procedure  is  taken  -for  all  elements.   Then  the  total 
stiffness  matrix  is  arrived  at  by  combining  the  element 
stiffness  matrices  using  structural  assembly  techniques 
(Zienkwicz,  1971) 

For  a  grid  o-f  elements,  a  system  of  banded  simultaneous 
equations  results: 

CS3   Ch>   =    CF>  (3.17) 

where  CS]  is  the  global  matrix  o-f  cof f icients  which 

incorporates  the  properties  o-f  the  materials  and  the  geometry  of 

the  elements;  Chy  is  the  vector  o-f  the  unknown  h's  at  the 

nodes;  and  CFy  is  the  load  vector.   After  the  introduction  of 

appropriate  boundary  conditions,  this  system  of  equations  can 

now  be  solved  by  computations  for  the  unknowns,  h,  yielding  the 

piezometric  head  distribution  in  the  solution  domain. 
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Free  Sur-face  Computation 

The  flow  domain  must  initially  be  established  assuming  an 
arbitrary  free  surface,  i.e.  the  user  initially  estimates  the 
free  surface  loacation.   All  boundary  conditions  are    solved  for 
the  preset  conditions  of  the  problem.   The  nodes  which  are    not 
fixed  are    allowed  to  move,  where  during  the  adjustment  the 
interior  nodes  are    also  moved  in  order  not  to  have  distortion 
of  the  elements.   Then  additional  iterations  are    made  until  all 
boundary  conditions  and  user  prescribed  error  correlations  have 
been  met  for  the  configuration. 
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3.3  Variables  and  Dimensional  Analysis 

The  variables  investigated  in  this  study  ar&    shown  in  Figure 
3.1.   For  the  analysis  o-f  seepage  using  the  -finite  element 
program,  variables  were  changed  -For  each  computer  run.   From 
each  run  data  was  generated  to  compare  the  changes  occur ing  -for 
each  variable  change. 

Only  the  two-dimensional  -flow  condition  was  considered  in  this 
study.   For  this  -flow  condition,  the  building  drainage  system 
is  in  one  cross-section  o-f  the  building  where  each  drain 
collects  the  seepage  and  collector  pipes  or  drains  take  the 
drainage  by  gravity  to  a  sump  or    some  exit  drain.   This  is 
depicted  by  Figure  3.4. 
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Figure   3.4     Underdrain  System  For   Seepage   Removal. 
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The  total  discharge  of  the  system,  qT  ,  can  be  express  as  the 
following  -Function: 

qT  =  j<    Hr,  h  ,  dw,  Sp,  K  ,  %9  8  )       (3.18) 
where      qT  =  total  discharge  of  all  drains,  i.e.  the  sum  of 

discharges  through  each  seepage  face  for  the 

model,  ft  J  /day/ft 
Hr   =  depth  of  saturated  soil  from  the  drain 

discharge  elevation  to  the  impervious 

boundary,  ft 
h    =  fixed  head  of  above  the  drains  ,  ft 
dw  =  drain  width,  ft 
Sp   =  drain  spacing,  ft 

K    =  hydraulic  conductivity  of  the  soil,  ft /day 
S  FS  =  difference  in  elevation  of  the  free  surface  to 

the  water  level  of  each  drain,  ft 
0   =  slope  of  the  impervious  boundary  with 

horizontal,  ft/ft 

Dimensional    analysis    of    equation     (3.18)     yields: 


q  T  f    I    H  d  S  .  (3. 19) 

Civ  w  p  \ 


K    h  *    h  h         h 


Since  permeability  is  moved  to  the  left  side  of  equation  (3.19) 

these  variable  is  made  proportional  to  the  discharge.   The 

elevation  of  the  free  surface,  SFS ,  i s  an  output  variable,  so 
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it  is  not  given  in  equation  (3.19.   Also  since  the  slope  of  the 
impervious  boundary,  9  ,  is  already  di mensi onl ess  it  can  be 
given  as   a  single  term  in  FT/FT. 

The  e-f-fects  of  these  functional  variables  and  other  results  are 
given  in  the  form  of  dimensi onless  plots  in  Chapter  4.   A 
discussion  of  the  results  is  given  in  that  chapter.   Since 
considerable  interest  is  in  quantities  of  discharge  and  the 
location  of  the  free  surface,  this  analysis  can  be  used  to  see 
trends  and  effects  that  each  variable  has  on  the  seepage 
domain. 
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3.4  Program  Veri f i cati on 

To  verify  that  the  computer  program  used  in  this  study  is 
giving  valid  results,  two  other  methods  of  solutions  were  used 
■for  a  selected  check  case.   Another  micro  computer  finite 
element  program,  SEEP,  and  an  appropriate  analytical  method 
based  on  Dupuit's  assumption,  was  used  for  the  simple 
configuration  of  a  drainage  blanket  with  a  single  pipe  type 
drain  (Duncan  and  Wong,  1985).   This  is  shown  in  Figure  3.5. 


/&s&y&//          \ 

.v  ■•■■..  ■.  •■  -.■..•■■-•  ■•■■.-  .-.•■  •■'..-•-.■*.-.  • 

Figure  3.5   Check  case  example  of  a  blanket  drain  with  single 
pipe  trench  drain. 


A  grid  was  generated  for  the  SEEP  program  to  match  the 
configuration  and  boundary  conditions  as  by  Figure  3.4.   The 
total  flow  for  this  check  case  is  calculated  by  the  SEEP 
program  as  based  on  the  Dupuit  assumption  of  a  fixed  domain  and 
fixed  amount  of  head.   Seep  also  tabulated  the  flow  quantities 

for  all  nodal  points  with  fixed  heads.   For  this  model,  the 
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■free  surface  converged  on  the  free  surface  location  with  an  error 
criteria  of  0.16  feet. 

The  location  of  the  free  surface  as  computed  by  the  Georgia 
Tech  model  converged  to  the  free  surface  having  an  error  factor 
of  only  0.001  feet,  a  much  more  accurate  and  sensitive 
convergence.   However,  the  two  models  compared  very  close  as 
far  as  the  location  of  the  free  surface.   Table  3.1  lists  the 
elevations  of  the  free  surface  nodes  for  the  two  models. 
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TABLE  3.1    Elevations  o-f  the  Free  Surface  Nodes 


x  coordinate 


SEEP:VPI 
Model 


Ga  Tech 
Model 


y  coordinate 


(I) 

40 

70 

80 

85 

87.5 

90 

92 

95 

97.5 
100 
106 
115 
1 30 
1  50 


78.00 
72.55 
67.07 
64.27 
62.03 
60.39 
58 .  00 
58.00 

60.  08 

61.  16 
61.92 
63 .  00 
63.76 
64.26 
64.42 


78 .  00 

72.73 

67.57 

65.  17 

63.84 

58.73 

58.  00^»-_  drain  location 

58.00 


58. 

62. 

15 

62. 

52 

63. 

23 

63. 

,93 

64. 

42 

64. 

62 
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To  veri-fy  that  -flow  quantities  were  accurate  -for  the  Ga  Tech 
model,  results  were  also  compared  with  the  results  obtained 
■from  a  Dupuit  assumption  solution,  developed  by  the  Federal 
Highway  Administration's,  Highway  Subdrainage  Design  Manual, 
(1980).   Figure  3.6  was  used  to  calculate  the  quantity  o-f  flow 
to  the  drain.   The  following  calculations  based  on  Figure  3„6 
are    presented: 


L  .     91' 

=  0.017   , =   -   1.57 


H      58' 
o 


H       o 
o 


erg. 


KCH  -  HQ)     0.2(78' -58' ) 


32 


92 


2.17   (-from  Figure  3.5) 


there-fore   q  2 


0.2(20) 


2.  17 


=   1.84   CFD 


q    1.84 
2_.  

K    0 .  2 


(H  -  H  )2 


—   7.  JLJL  , 


_Q_ 


K     2(L  -b) 


(20) 


(91-1) 


2  22 


q  d 

—  =    9.22  +  2.22    =    11.44 

K 

ther-fore    qd  -  0.2(11.44)  -   2.29    CFD 

The  results  o-f  this  computation  is  compared  with  the  results 
o-f  the  two  -finite  elements  in  Table  3.2 
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Figure  3.6    Chart  for  Determining  Flow  Rate  in  Symmetrical 
Underdrains   (after  FHA,  1980) 
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Table  3.2   Discharge  to  the  drain  check  cases,  CFD 

Discharge,  CFD 
SEEP: VPI  Model  2. 15 

Ga  Tech  Model  2.34 

FHA  Nomograph  2-29 

In  summary,  the  Ga  Tech  model  compares  very  well  with  other 
model  and  with  the  FHA  nomograph  -for  -flow  quantities.   The  Ga 
Tech  model  was  -found  to  more  error  criteria  sensitive.   The 
•free  surface  location  check  also  proved  to  compare  very  well 
with  the  SEEP  model. 
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CHAPTER  4 

RESULTS  AND  DISCUSSION 

4.1  Introduction 

This  chapter  discusses  the  analysis  and  results  of  the  two 
dimensional,  steady  state,  free  surface  finite  element  study  of 
seepage  to  building  underdrain  systems.   Details  presented 
include   model  set-up  such  as  boudary  conditions,  grid 
development  and  basic  first  assumptions.   Output  variables 
considered  are    total  discharge  and  discharge  at  each  drain, 
plus  the  free  surface  elevation  between  the  drains. 

In  studying  sensitivity  of  the  model  for  the  input  variables, 
it  was  found  that  the  distance  from  the  Dirichlet  boundary 
face,  or  constant  head  boundary,  Li,  to  the  first  drain  was 
critical.   Li  was  found  to  be  critical  in  terms  of  output 
charcteristics,  i.e.  discharge  and  location  of  the  free  surface 
between  drains.   Therfore,  two  sets  of  results  are    presented  for 
two  Li  distances  of  400  feet  and  80  feet. 

4.2  Model  Set-up 

Figure  4.1  shows  the  basic  grid  that  was  developed  for  this 

seepage  problem.   Using  150  elements  and  521  nodes,  the 

geometry  of  the  mesh  was  based  on  a  five  drain  system.   The 
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sides  and  bottoms  of  the  drains  were  modeled  as  seepage  -Faces. 
Using  the  two-dimensional  model,  the  fifth  drain  is  considered 
to  be  at  the  building  centerline.   Assuming  symmetry,  only 
one-half  of  the  building  was  analysed.   The  left  most 
boundary,  as  shown  in  Figure  4-1,  is  a  Dirichlet  boundary 
condition  which  has  the  constant  head  required  to  allow  for  the 
seepage  to  enter  the  drains. 

Listed  below  are    some  basic  assumptions  for  the  study  made  for 
this  seepage  problem: 

<1)   Steady  state  flow  conditions 

(2)  Two-dimensional  analysis  where  the  right  most 
boundary  is  the  building  centerline  and  by  symmetry  the  flows 
are  only  one-half  of  the  total  for  the  total  width  of  the 

bui lding. 

(3)  Free  surface  nodes  move  to  the  computed  phreatic 
level,  i.e.  piezometric  head  equals  atmospheric  pressure  at  the 
upper  most  nodes 

(4)  Partial  saturated  flow  is  not  considered 

(5)  Evaporation  is  not  considered 
<6)  Retardation  is  not  considered 
(7)   Infiltration  is  not  considered 

(3)   The  soil  is  homogeous  being  either  isotropic  or 
ani  sotropic. 
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4.3  Variables  Used 

For  this  study,  key  variables  were  used  in  the  -finite  element 
model  to  look  at  their  effects  on  the  -flow  of  seepage  -for  the 
building  underdrain  system.   Two  major  design  parameters 
necessary  to  model  any  seepage  problem  is  the  total  discharge 
of  the  drain  system  and  the  location  of  the  free  sur-face.   For 
this  study,  the  -following  output  variables  ^re    used: 

qd   =   discharge  or  seepage  quantity  to  each  drain 

(cubic  -feet  per  day  per  linear  foot) 
q^   =   total  discharge  of  all  drains  (same  units) 
S  pC  =   height  of  free  surface  above  drain  elevation. 

Discharge  quantities  are  in  units  of  flow  quantity  in  time  per 

unit  length  of  the  total  flow  continum.   Therfore,  flow  for  a 

total  system  should  be  multiplied  be  the  length  of  the 
structure,  i.e.  in  the  third  dimension. 

Input  variables  used  in  this  problem  to  study  their  effects  on 
seepage  are    listed  below: 

Li   =   distance  from  the  left  most  constant  head 

boundary  to  the  first  drain  (feet) 

K    =   hydraulic  conductivity  (feet /day) 

Kv   =   hydraulic  conductivity  in  the  vertical 
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direction  (feet /day) 
Kh   =   hydraulic  conductivity  in  the  horizontal 

direction  (-feet  /day) 
h    =   hydraulic  head  of  -flow  continum  above  the 

drain  elevation  (feet) 
Hr   =   depth  of  aquifer  between  the  elevation  of  the 

drain  and  the  impervious  boundary  (feet.) 
Sp   =   spacing   between  drains  (feet) 
Dw   =   drain  width  (feet) 
6    =   slope  of  the  impervious  boundary  (7. ) 


For  better  reference  to  the  notation  of  the  variables  refer  to 
Figure  3.1  in  Chapter  3. 

For  this  f  i  n  i  t  e  e  1  e m e n  t.  p r  o b  1  e m  t  h e  u. nits  of  fee t  a n d  cl a y s  w e r  6 
used  to  give  results  for  the  output  variables  in  cubic  feet  per 
day  and  feet  per  day  for  velocity,, 

H y draulic  c  o n d  u c  t  i v  i  t  y  o  r  p e r  m e  a b  i 1  it  y  v  a 1 u  e  s  u  s e d  f o r  t  h x  s 
study  were  0.0002,  0.02,  0.2,  and  4  feet  per  day,,   A  soil's 
hydraulic  conductivity  having  a  value  of  0.0002  feet  per  day  is 
equivalent  to  a  fine  grained  soil  such  as  a  clay.   Hydraulic 
conductivities  are)    tabulated  in  Table  4.1  to  show  cross 
reference  to  soil  types  and  permeabilities  in  units  of 
c  en  t  i  meter  s  p er  sec  on  d . 
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T a b  1  e  4 .  .1.  H y cl r anlic  C a n d u c t  i  v i  t y  V a  1  ues  F a r  V a rio u s  So i  .1  T y pes 


f-ine  g r  a  i  ned  so i  I 
Fine  grained  soil 
silt 

sand , si  1 1 
clean  sand 


0.0002  feet/dav 


0.002  feet /day 


0.02  feet/dav 


0.2  feet/day 


-7 
/:•: 1U   cm /sec 


7>:  10   cm /sec 

7:-;  10   cm /sec 

-7 
7x10   cm /sec 

-7 
4  -feet/dav     1.4x10   cm/ sec 


The  other  input  variables  used  were  typically  realistic  values 
to  meet  typical  field  conditions.   Listed  below  are    some  of 
this  values,  where  for  each  computor  run  specific  variables 
were  changed  and  the  others  held  constant.   The  intent  was  to 
see  what  effects  that  each  variable  change  made  on  the  model. 

Li  =  80' ,  400'  and  1000' 

h   =  10' ,  20'  and  40' 

Hr  =  20' ,  5B'  and  100' 

Sp  -  15' ,  25'  and  35' 

Dw  =  1 ' ,  2'  and  4' 

9   =  0"/.,  57. ,  107  and  207. 

Kv  =  .0002,  0.2,  4  and  15 

Kh  =  0.2 

K   =  0.0002,  0.02,  4 
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4.4   Model  Sensitivity 

The  model  was  developed  assuming  a  homogeneous  isotropic  media, 
except  -For  a  -few  cases  where  anisotropic  conditions  were  studied. 
The  nodes  of  the  top  three  rows  o-f  elements  were  allowed  to  move 
■for  iteration  o-f   the  free  surface.   In  all  runs,  the  maximum 
number  of  iterations  was  six  and  on  an  average  three  to  four  were 
found  to  be  adequate.   The  model  was  made  to  iterate  to  an 
accuracy  of  0.001  feet  during  the  free  surface  moves. 

In  calibration  of  free  surface  models  careful  attention  must  be 
made  to  the  moving  nodes.   Problems  exist  when  the  elements  become 
less  than  zero  in  area  or  when  the  elements  actual  may  become  so 
distorted  that  they  may  become  inverted.   This  was  not  the  cas& 
for  the  mesh  used  in  this  study. 

In  studying  the  flow  quantities  and  velocities  for  the  basic 

geometries  and  variable  changes,  all  results  seemed  reasonable  and 

observed  trends  appeared  physically  possible.   However  large 

variation  of  output  was  noticed  by  the  change  in  distance  from  the 

constant  head  boundary  to  the  first  drain,  Li.   The  first  set  of 

results  is  presented  based  on  the  influence  length,  Li,  being  10 

times  the  maximum  head  difference,  i.e.  10  x  40'  =  400'.   Depth  of 

aquifer  for  this  case  was  held  constant  at  100  feet.   Use  of  10 

times  the  head  differential  is  based  on  Bear's  (1979) 

recommendat i on . 
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Figure  4.2  depicts  the  moved  -free  surface  corner  nodes  for  an 
influence  length  of  400  feet.   For  comparison.  Figure  4.3  shows  thi 
moved  free  surface  corner  nodes  for  Li  equal  to  80  feet  and  a 
thickness  of  aquifer,  Hr ,  equal  to  58  feet. 
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Figure  4.4  and  4-5  are    presented  to  show  the  dramatic  influence 
that,  the  distance  -from  the  constant  head  boundary  from  the  -first 
drain,  Li,  has  on  total  discharge  and  the  maximum  -free  surface 
elevation,  respectively.   For  further  reference,  the  total 
discharge,  q  ,  is  the  sum  of  the  drain  discharges  and  the  maximum 
free  surface  is  the  highest  elevation  of  water  above  the  drain 
elevation.   Lengths  between  300'  and  400'  appear  to  result  in  less 
change  in  flow  quantity  and  free  surface  changes.   This  equates  to 
the  distance  where  the  constant  head  boundary  away  from  the  drains 
begins  to  cause  the  least  effect  on  discharge  and  free  surface 
changes  is  approximately  15  to.  20  times  the  head  differential. 
This  perhaps  is  good  evidence  that  the  actual  field  conditons  for 
each  particular  application  of  groundwater  movement  should 
carefully  be  looked  at  in  terms  of  this  distance. 
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Figure  4.4   Distance  -from  the  constant  head  boundary  to  the  -first 
drain  vs.  total  discharge  of  drains, 
h  =  20',  Sp  =  IS',  K  =  0.2  -ft /day 
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Figure    4.5    Distance    -from    the    constant    head    boundary    to    the    first 
drain    vs.     maximum    -free    surface    elevation, 
h    =    20',     Sp    =    IS7,     K    =    0.2    -ft/day 
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4.5  Results 

This  section  will  present  the  computer  results  in  the  -form  of 
di  mensi  onl  ess  plots  and  other  nomographs  to  show  the  effects 
that  each  variable  change  has  on  the  seepage  problem.   Two 
sets  o-f  plots  are    given  where  the  first  set  is  -for  the  constant 
head  boundary  distance,  Li,  equals  400  -feet  and  the  depth  of 
rock,  Hr ,  equals   100  feet.   The  second  set  of  plots  are 
essentially  the  same  as  the  first  except  Li  equals  80  feet  and 
the  Hr  equals  58  feet.   For  variable  notation  refer  to  Figure 
3.  1. 

Drain  Discharge  (Li  =  400J  and  Hr  =  100') 

Figure  4.6  shows  the  influence  of  drain  width,  Dw,  for  head 
differentials,  h,  of  10,  20  and  40  feet  on  total  discharge  of 
the  drains.   This  shows  that  drain  width  is  not  particularly 
critical   for  low  head  differences.   It  is  reasonable  to  assume 
that  total  discharge  could  approach  a  constant  for  drain  widths 
greater  that  four  feet  and  head  less  than  20  feet.   Although 
total  discharge  increases  with  drain  width,  this  plot  does  not 
show  what  drain  width  is  the  most  effective  or  most  efficient. 
However,  the  influence  of  head  can  be  easily  observed 
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Dw 


Figure  4„£>   Total  drain  discharge  vs.  drain  width  with  head 

di  f -f  erenti  al  s. 

Li  =  400*,  Hr  =  100*,  Sp  =  25*,  K  =  0.2  ft /day 
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Figure  4.7  depicts  the  results  a-f  varying  anisotropic:  soils 
with  varying  heads.   In  this  analysis,  vertical  hydraulic 
conductivity   varies  while  horizontal  conductivity  was  held 
constant.   Figure  4.7  shows  that  with  increase  in  vertical 
hydraulic  conductivity  discharge  increases.   This  appears  to 
satisfy  theory  ,  that  with  less  resistance     -flow  increases. 
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Figure  4.7   Total  drain  discharge  vs, 
head  di-f-f erential . 
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Figures  4.8  through  4.12  are    graphs  which  show  the  discharge 
quantities  -For  each  o-f  the  -five  drains.   Figure  4.8  depicts  the 
effects  o-f  depth  o-f  impervious  layer,  Hi- ,  on  the  drain 
discharge  -for  each  o-f  the  drains.   This  data  shows  that  with 
increased  depth  discharge  increases.   It  is  interesting  to  note 
that  practically  no  -flow  is  recorded  for  drain  5  when  a  shallow 
depth  of  impervious  layer  exists,  i.e.  Hr  =  20  feet. 
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Figure  4.8   Drain  discharge  vs.  thickness  of  aquifer. 
Li  =  400',  h  =  20',  Sp  =  25 ? ,  Dw  -  27  ,  K  ■  0.2  ft/day 
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Figure  4.9  shows  the  effects  of  head  di  f  f  erent  i  al  on  the 
discharge  of  each  of  the  drains.   It  is  interesting  to  see  that 
though  the  first  drain  is  the  most  critical  drain,  i .  e.  a 
large  amount  of  the  total  flow  is  collected  in  this  drain,  the 
remaining  drains  collect  discharge  decreasing  linearly  from 
the  left.   This  is  possiblly  due  to  the  influence  of  the 
constant  head  boundary  distance,  Li,  being  400  feet  away.   This 
shows  that  for  the  drains,  flow  quantity  increases  linearly  with 
increase  in  head. 
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Figure    4.9      Drain    discharge    vs.     head    di f -f erenti al  . 

Li     =    400',     Hr    =    100*,     Sp    =    25',     Dw    =    2',     K    -    0.2    -Ft/day 
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Figure  4,10  depicts  the  results  of  comparing  the  drain 
discharge  versus  the  spacing  between  drains.   It  is  of 
interest  to  see  that  drain  number  one  -For  spacing  equal  to  25 
■Feet,  collects  more  seepage  than  the  other  discharges  for  the 
other  spacings.   However,  at  drains  3  through  5  this  changes 
and  the  discharge  for  15  feet  spacing  is  greatest.   This  will 
be  analyzed  later  in  the  dimensional  analysis. 
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Figure  4.10   Drain  discharge  vs.  drain  spacing. 
Li  =  400',  Hr  =  100',  h  =  20%  Dw  =  2?  ,  K  =  0.2  ft/day 
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Figure  4.11  shows  the  e-f-fects  of  ani  sotropy  of  the  soil 
matrix  on  drain  discharge-   It  is  of  interest  to  see  that 
the  discharge  -for  Kv/Kh  equal  one  is  much  greater  than  the 
other  drain  discharges  -for  Kv/Kh  less  than  one. 
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Figure  4.11   Drain  discharge  vs.  anisotropy. 
Li  -  400*,  Hr  =  100%  Sp  =  15',  Dw  =  2,  h  =  20' 
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Figure  4.12  depicts  data  generated  -for  the  condition  o-f  a 

sloping  impervious  layer.   The  slope  o-f  the  boundary  is  given 

as  a  percentage  o-f  Ax/ty.        The  depth  o-f  the  layer  is  held 

constant  at  drain  one  and  the  layer  slopes  clockwise  -from  the 
positive  x  horizontal  axis. 

Though  actual  quantities  o-f  discharge  are    only  changing  in  very 
small  ammounts,  the  trend  is  seen  that  -for  slopes  greater  than 
5"/.  the  discharge  increases  in  drains  3  through  5.   However  at 
some  slope  between  57.  and  10"/.  the  discharge  decreases.   It  is 
also  seen  that  the  discharge  of  drain  -five  is  the  least  e-f-fected 
by  any  slope   change,  in  terms  o-f  discharge  quantity. 
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Figure  4.12   Drain  discharge  vs.  slope  o-f  the  impervious 
1 ayer - 

Li  =  400*,  Hr  =  100'  (-first  drain),  h  =  20 ?  ,  Sp  =  25",  Dw 
K  =  0.2  -ft /day 


76 


Free  Surface  Location  (Li  =  400'  and  Hr  =  100') 


Figure  4.13  shows  the  in-fluence  of  the  drain  width  on  the 

maximum  -free  surface  height,  S  ttc    u   In  all  cases  this  occured 

*   •    "max 

between  drains  one  and  two.   This  plot  shows  how  the 
relationship  is  linear  and  that  the  height  o-f  the  -free  surface 
is  not  greatly  affected  by  the  drain  width. 
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Figure  4.13   Maximum  -free  sur-face  height  vs.  drain  width  with 

varying  heads. 

Li  =  400*,  Hr  =  100%  Sp  =  257  ,  K  =  0.2  ft/day 
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Figure  4.14  depicts  the  influence  of  drain  spacing  and  head  on 
the  maximum  -free  surface  height.   These  results  show  that 
increased  drain  spacing  causes  an  increase  in  the   free 
surface.  Therefore  one  should  consider  the  smaller  spacing  or 
possibly  deeper  drains. 
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Figure  4.14  Maximum  free  surface  height  vs.  spacing  of  drains 

with  varying  head. 

Li  =  400',  Hr  =  100',  Dw  =  2' ,  K  =  0.2  ft/day 
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Figure  4.15  shows  the  relationship  of  maximum  -free  surface 
heioht,  Sv<;   *  as  a  faction  of  head,  h,  and  anisotropy,  Kv/Kh 
The  same  trend  is  seen  here  as  noticed  with  discharge  for 
anisotropy,  that  with  increased  vertical  hydraulic 
conductivity  the  height  of  the  free  surface  increases. 
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Figure  4.15   Maximum  free  surface  height  vs.  anisotropy  with 

varying  head. 

Li  =  400,  Hr  -  100*,  Sp  =  25 '  ,  Dw  =  2' 
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In  Figures  4.16  through  4.20,  the  free  surface  corner  nodes  are 
plotted  -for  the  different  cases  of  analysis.   To  see  the 
effects  of  the  free  surface  changes  versus  the  variable  changes 
only  nodes  between  the  drains  are    plotted.   In  the  sence  of 
design,  the  max i mum  height  of  the  free  surface  would  be  of  the 
most  concern,  which  is  in  between  the  first  two  drains. 
However,  it  is  of  interest  to  see  the  total  free  surface  from  a 
standpoint  of  total  cross  section  of  the  building. 

In  Figure  4.16  the  effect  of  thickness  of  aquifer  on  location 

of  the  free  surface  is  shown.   It  is  of  interest  to  see,  as  in 

the  case  for  discharge  shown  by  Figure  4.8,  that  for  shallow  depths 

of  impervious  layer,  Hr  =  20 * ,  that  very    little  rise  in  the 

free  surface  occurs  in  between  drains  four  and  five. 

Figure  4.17  shows  the  effects  of  head  differential  on  the 
location  of  the  free  surface  nodes.   It  is  of  interest  to  see 
for  cases  of  head  equal  to  20  and  10  feet  the  free  surface 
nodes  are  approaching  gentle  mounds  of  water  in  between  drains. 
However  for  the  head  equal  to  40  feet  the  free  surface  peaks 
rather  abruptly  between  the  first  two  drains  but  resumes  the 
gentle  fall  in  free  surface  between  the  other  drains. 

Figure  4.18  shows  the  effects  of  drain  spacing  on  the  location 

of  the  free  surface  nodes.   This  shows  as  spacing  is  increased 

the  free  surface  increases. 
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For  Figure  4, 19,  the  effects  of  the  sloping  aquifer  on  location 
of  the  free  surface  nodes  is  shown.   The  same  trend  is  seen 
here  as  in  discharge,  shown  by  Figure  4.12,  that  for  slopes 
greater  than  10'/.  the  free  surface  location  drops.   The  change 
in  free  surface  nodes  is  very  small  for  the  cases  of  slope 
equal  to  5"/.  and  10"/.. 
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Dimensionsal  Analysis  (Li  =  400'  and  Hr  =  IPO"  ) 

Figures  4.21  through  4.24  are    plots  o-f  the  data  generated  by  a 
dimensional  analysis  as  discussed  in  Chapter  3.   The  discharge 
values  given  are    unit  discharge  quantities  and  each  case  is  -for 
the  particular  variables  noted  in  each  -figure.   Unit  discharge 
is  de-fined  as  the  sum  o-f  the  total  drain  discharges  for  each 
o-f  the  drains  per  unit  width  in  the  z  axis  or  the  third 
di  men si  on . 

In  Figure  4.21  the  effect  of  drain  width  is  analyzed.   This 
plot  shows  that  with  increased  drain  width  discharge  increases, 
but  with  increased  head  the  discharge  may  increase  at  a  much 
faster  rate.   From  this  analysis  it  would  appear  that  the 
greatest  increases  of  flow  rate  would  be  for  the  drain  widths 
between  one  and  two  feet  wide.   The  following  analysis  is 
performed: 


Drain  Widths 


head 


10  0.034  0.053  0.065  0.073 
increase  —  0 .019     0 . 0 12       0 . 009 

20  0.035  0.055  0.070  0.086 
increase  »-  0.02      0.015       0.016 

40     0.036     0.057      0.075       0.09 

increase  —  0.021     0.018       0.015 
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This  analysis  shows  the  increases  o-f  the  -flow  -function  -for  the 
di-f-ferent  head  increases  as  shown.   The  maximum  increase  is 
seen  for  the  values  between  drain  widths  o-f  one  and  two. 
Ther-fore  the  most  e-f-fective  increase  in  -flow  occurs  -for  drain 
widths  between  1  and  2  -feet.   Another  relationship  could  be 
stated,  that  -for  Dw/h  ratios  less  than  0.2  the  -flow  quantity 
increases  at  a  rapid  rate. 
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Figure    4.21       Dimensional    analysis,     Drain    width, 
Li    =    400',     Hr    =    100',     Sp    =    25',    K    =    0.2    -ft/day 
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In  Figure  4.22  the  effects  o-f  thickness  of  aquifer  is  shown  -for 
the  dimensional  analysis.   This  plot  concludes  that  -for  Hr/h 
ratios  greater  than  2  or  3  the  net  increase  in  -flow  is 
approaching  a  linear  relationship.   This  states  that  with 
increasing  thickness  o-f  the  aquifer  the  dicharge  is  expected  to 
increase  very  near  a  linear  increase  and  for  very  shallow 
thicknesses  the  flow  becomes  very  small. 
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Figure  4.22   Dimensional  Analysis,  Depths  of  impervious  layer 
Li  =  400',  Sp  =  25?,  Dw  =  2',  K  =  0.2  ft/day 
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The  effects  of  drain  spacing  is  presented  in  Figure  4.23.   The 
■flow  function  -for  this  analysis  is  treated  di  f  f  erent  1  y  -For  the 
spacing  changes.   As  the  spacing  was  increased  -for  each  o-f  the 
computer  runs  the  flow  domain  increased  accordingly  to 
accommodate  the  five  drains  at  the  appropriate  spacing. 
Therfore  to  gain  understanding  of  the  effects  of  flow 
relationships  with  respect  to  drain  spacing,  each  term  was 
normalized  by  the  corresponding  spacing  between  the  drains  for 
each  head  increase.   This  analysis  shows  for  decreased  spacing 
discharge  increases  and  as  head  increases  so  the  discharge. 
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Figure  4.23   Dimensional  Analysis,  Drain  spacing. 
Li  =  400,  Hr  ■  100'  ,  Dw  -  2'  ,     K  -  0.2  -ft /day 

Figure  4.24  shows  the  dimensional  analysis  -for  ef-fects  of  the 
sloping  aqui-fer.   For  this  particular  -flow  matrix  as  modeled, 
a  de-finate  break  in  -flow  increase  is  noticed.   The  peak 
discharge  -function  is  at  slope,  0   equal  0.05. 
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Figure  4.24  Dimensional  Analysis,  Slope  of  aquifer. 
Li  =  400' ,  Hr  =  100'  (-first  drain),  Sp  =  25'  ,  Dw  =  2' 
■ft/day 


=    0.2 
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Results  For  Li  =  80'  and  Hr  =  58' 

The  following  plots  and  nomographs  are    the  results  for  the 
computor  runs  where  distance  from  the  constant  head  boundary  to 
the  first  drain,  Li,  equals  BO  feet  and  the  thickness  of  the 
aquifer  equals  53  feet.   The  plots  are  essentially  the  same  as 
for  the  previous  case  except  these  variables  are    held 
constant . 

Drain  Discharge  (  Li  =  BO'  and  Hr  =  58'  )  i 

Figure  4.25  shows  the  effects  of  drain  width  on  total  drain 
discharge.   The  results  are    generally  the  same  as  shown  in 
Figure  4.6  except  discharge  quantities  range  from  3  to  4  times 
greater. 
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Figure  4.25   Total  drain  discharge  vs.  drain  width  with  varying 

head. 

Li  =  SO',  Hr  =  58*,  Sp  =  25*,  K  =  0.2  ft/day 
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The  results  of  anisotropy  of  the  soil  continum  is  shown  by- 
Figure  4.26  for  the  short  Li  distance  of  80  feet.   The  trend  is 
slightly  different  than  for  the  case  for  Li  equal  400  feet  as 
shown  by  Figure  4.7.   For  the  largest  head  value  of  40  feet  the 
discharge  decreases  just  slightly  with  decrease  in  vertical 
hydraulic  conductivity,  but  for  a  Kv/Kh  ratio  less  than  0.1  the 
discharge  decrease  is  very  dramatic.   However  for  the  low  head 
value  of  10  feet  the  general  decreasing  trend  is  prevelant  as 
the  case  for  Figure  4.7. 

For  this  case  the  flow  quantities  have  increased  approximately 
4  to  5  times.   This  shows  the  large  effects  that  the  distance 
of  the  constant  head  boundary  has  on  the  seepage  quantities. 
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Figures  4.27  and  4.28  depict  the  same  general  trend  o-f  maximum 
^low  at  drain  one  and  a  general  decrease  in  -flow  tor  each 
consectutive  drain  is  noticed  -for  e-f-fects  o-f  thickness  o-f 
aqui-fier  and  head  respectively.   Flow  quantities  are    again  3  to 
4  times  greater  -for  this  case  than  -for  the  earlier  case  as 
shown  by  Figures  4.8  and  4.9.   For  this  case  of  shorter 
distance  to  the  constant  head  boundary  equal  80  -feet,  the 
mid-peak  discharges  at  the  third  drain  is  not  noticed  as  in 
Figure  4.9.   In  Figure  4.28  it  is  oi     interest  to  see  that  the 
largest  head  value  equal  40  -feet  has  such  a  large  e-f-fect  on  the 
•first  drain  collecting  a  large  portion  of  the  total  discharge 
o-f  the  system. 
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Figure  4.27   Drain  discharge  vs.  depth  o-f  aqui-fer. 
Li  =  80',  h  =  20',  Sp  =  25',  K  =  0.2  -ft /day,  Dw  ■ 
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Figure  4.28   Drain  discharge  vs.  head  differential 
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The  effects  of  drain  spacing  for  distance  of  the  constant  head 
boundary  to  the  first  drain,  Li,  equal  80  feet,  is  shown  by 
Figure  4.29.   The  same  trends  are    seen  for  this  case  as  that 
shown  by  Figure  4.10.   Flow  increase  is  somewhat  less  in  this 
case  compared  to  the  earlier  case  of  Li  equal  400  feet.   Flow 
quantities  appear  to  range  from  2  to  3  times  as  great. 
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Figure  4.29   Drain  discharge  vs, 
Li  -  SO7,  Hr  =  58'  ,  h  =  20',  Dw 
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In  Figure  4.30  the  effects  of  slope  of  the  impervious  boundary, 
°  ,  is  shown  -for  drain  discharge.   A  significant  change  is 
noticed  for  this  plot  when  compared  to  the  earlier  plot  in 
Figure  4.12  for  sloping  aquifer.   For  Li  equal  to  400  feet,  the 
peak  discharge  occured  for  a  slope  between  5  and  10  percent. 
In  this  case,  where  Li  equals  80  feet,  the  peak  flow  is  for  the 
case  of  the  steepest  slope  equal  to  20  percent.   For  drain  one, 
the  flow  is  8  times  greater  than  the  flow  for  drain  one  as 
shown  by  Figure  4.12,  Li  -  400  feet.   The  flow  is  5  times 
greater  in  drain  five  in  this  case  than  the  flow  in  drain  five 
for  the  case  of  Li  =  400  feet. 
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Figure  4.30   Drain  discharge  vs.  slope  of  the  impervious 

boundary. 

Li  =  80',  Hr  =  SB'  (-First  drain),  h  =  20 '  ,  Sp  =  25'  ,  Dw  =  2? 

K  =  0.2  ft /day 
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Free  Surface  Location  (  Li  =  SO*  and  Hr  =  5Q-      ) 

Figures  4.31  through  4.36  show  the  results  of  the  -free  surface 
location  plots  -for  the  case  o-f  Li  equal  80  -feet  and  Hr  equal  58 
•feet.   The  same  general  trends  ar&    shown  for  this  plots  as 
shown  by  the  earlier  case  of  Li  equal  400  feet  and  Hr  equal  100 
feet.   In  general,  the  free  surface  has  moved  up  approximately 
three  times  in  height.   An  exception  to  the  earlier  case  of  is 
shown  in  Figure  4.36.   This  figure  shows  the  location  of  the 
free  surface  nodes  with  the  effects  of  slope  of  the  impervious 
boundary.   It  is  somewhat  different  in  this  case  than  in  the 
earlier  case  as  shown  in  Figure  4.19  for  Li  equal  400  feet.   In 
this  case  the  free  surface  continues  to  increase  in  height  as 
the  slope  of  the  impervious  boundary  increases. 
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Figure    4.31       Maximum    -free    surface    height    vs.     drain    width 
Li     =    80',     Hr    =    58'.,     Sp    =    25',     K    =    0.2    ft/day 


105 


j .  b 


28   25   38 
HEAD  (FT) 


SPACIN 


VI 


1 1* 


SPACING 
251 


SPACING 
35' 
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Dimensional  Analysis  (  Li  =  SO'  and  Hr  =  58*  ) 

The  effects  of  drain  width  as  shown  by  the  dimensional  analysis 
is  shown  by  Figure  4.37.   The  general  trend  of  the  discharge 
increasing  with  drain  width  is  seen  in  the  data.   When  compared 
to  the  earlier  case  -for  Li  equal  to  400  -feet,  it  is  seen  that 
total  discharge  is  greatly  effected,  by  a  factor  of  three 
ti  mes.    H  ,  - 
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Figure  4.37   Dimensional  Analysis,  drain  width 
Li  =  80?,  Hr  =  58',  Bp  =  25',  K  =0.2  ft/day 
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Figure  4.38  shows  the  dimensional  analysis  results  -for 
thichness  of  the  aquifer.   In  this  case,  Li  equal  80  feet,  the 
relationship  is  more  linear  than  in  the  case  of  Li  equal  400 
feet  as  shown  in  Figure  4.22.   This  results  in  the  conclusion 
that  the  effects  of  Hr  has  a  pronounce  influence  on  the  flow 
quantities  when  coupled  with  the  distance  Li.   For  a  shorter 
distance  of  Li,  the  effects  of  depths  of  Hr  is  more  crictical 
and  more  of  a  linear  relationship.   However  when  the  distance 
from  the  constant  head  boundary  is  larger,  as  i  n  the  earlier- 
case,  the  effects  of  Hr  are  not  as  effective  when  Hr  i s  very 
smal 1 . 
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Fi  gure  4. 38 

Li  =  SO' ,  Sp 


Dimensional  Analysis,  depth  of  aquifer 
-25',  Dw  =  2'  ,  K  =  0.2  ft /day 
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Figure  4.39  shows  the  results  o-F  the  dimensional  analysis  -for 
drain  spacing.   Again  the  ordinate  is  normalised  with  respect 
to  the  area  subtended  by  each  -flow  continum  -For  each  case  o-F 
spacing.   Flow  quantity  would  equate  to  an  increase  o-F 
approximately  5  to  6    times  -from  this  case  of  Li  equal  80  -Feet 
versus  Li  equal  400  as  shown  in  Figure  4.23. 


Hr/h 


Figure  4.39   Dimensional  Analysis,  drain  spacing 
Li  =  BO',  Hr    =  5S7 ,  Dw  =  2? ,  K  =  0.2  ft/day 
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Figure  4.40  shows  the  dimensional  analysis  results  for  the 
sloping  impervious  boundary.   For  this  case,  Li  equal  80  -feet, 
as  depicted  in  the  earlier  plots  -for  sloping  aquifer  -for  drain 
dischrge  and  -free  surface  location,  the  trend  is  seen  that  for 
and  increase  in  slope  the  discharge  increases. 
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6 

Figure  4.40   Dimesional  Analysis,  sloping  impervious  boundary 
Li  =  80',  Hr  =  58' ,  Sp  -    25',  Dw  =  2,  K  =  0.2  ft/day 
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4.6   Summary 

This  chapter  has  included  the  results  of  many  computer  runs  -for 
the  analysis  of  a  building  underdrain  system.   The  importance 
of  model  sensitivity  must  be  stressed,  in  that  key  variables  may 
effect  the  results  very  significantly.   In  this  model,  it  was 
found  that  the  most  crictical  variables  affecting  the  total 
discharges  of  the  drains  and  location  of  the  free  surface  was 
the  distance  from  the  constant  head  boundary  to  the  drain 
system,  Li,  and  thickness  of  the  aquifer,  Hr . 

The  results  obtained  appear  to  be  reasonable  and  within  limits 
of  possible  field  occur ances.   However,  using  these  results  for 
quantifying  flows  and  the  free  surface  locations  requires  a 
very  careful  study  and  background  investigation  of  all 
geometric  configurations  and  soil  conditions.   It  is  intended 
that  this  study  has  presented  the  results  in  a  manner  that  one 
can  observe  the  effects  of  the  variables  studied  on  such  a 
underdrain  system. 

In  summary,  the  following  observations  are    discussed  for  the 
variables  used  in  this  analysis. 

Constant  Head  Boundary  Distance,  Li 


As  distance  from  the  constant  head  boundary  to  the  system 
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increases  the  flow  quantities  increase  -for  the  drain  system. 
In  terms  o-f  a  correlation  with  hydraulic  head  differential 
above  the  drain,  as  used  by  other  authors,  this  model  study 
shows  that  for  a  distance  of  approximately  15  times  the  head, 
discharge  appears  to  decrease.   In  terms  of  gradient,  this 
equates  to  a  1  over  15  ratio  or  0.067  gradient.   A  slightly 
higher  distance  was  observed  in  the  analysis  of  this  influence 
to  location  of  the  free  surface. 

Depth  of  Aquifer,  Hr 

The  results  show  that  with  increased  aquifer  thickness  the 
drain  discharge  increases,  increasing  very  near  to  a  linear 
relationship.   It  was  also  observed  that  for  the  very  shallow 
thicknesses,  the  last  drains  collected  very  little  flow.   The 
height  of  the  free  surface  is  affected  in  the  same  fashion. 
With  increased  thickness  of  aquifer  the  height  of  the  free 
surface  increases  between  drains. 

Drain  Spacing,  Sp 

Drain  spacing  effects  on  discharge  was  greatly  affected  by  the 

distance  of  the  system  from  the  constant  head  boundary.   A 

substantial  increase  of  5  to  6  times  the  discharge  was  observed 

for  the  two  distances  used  in  the  results  presented.   It  was 

observed  that  as  head  increases,  discharge  increases  for  all 

drain  spacings;  and  as  drain  spacing  increases,  discharge 
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decreases.   It  was  also  observed  that  -from  the  dimensional 
analysis,  the  most  efficient  drain  spacing  is  -for  a  spacing 
over  head  ratio  of  less  than  2.   The  location  of  the  free  surface 
was  observed  to  increase  between  drains  as  drain  spacing  is 
increased. 

Drain  Width,  Dw 

With  increased  drain  width  the  discharge  increases  for  all 
other  variables  constant.   However,  there  must  be  a  trade  off 
for  efficiency  of  the  complete  system.   It  appears  that  a  drain 
width  between  1  to  2  feet  is  the  most  effective  width.   This  is 
also  a  plus  from  a  standpoint  that  lateral  seepage  can  be 
easily  drained  into  thinner  trenches  or  deeper  constructed 
drai  ns. 

Slope  of  Impervious  Boundary,  0 

For  the  distance  of  the  constant  head  boundary  to  the  drain 

system  equal  to  80  feet,  the  drain  discharge  was  observed  to 

increase  with  increased  slope.   However,  for  the  distance,  Li 

equal  to  400  feet,  the  maximum  discharge  occurs  at  a  slope 

between  5  to  10  percent.   This  occur ance  is  not  fully 

understood  and  perhaps  should  be  further  studied  to  develope 

some  reasoning.   This  requires  that  careful  consideration  of 

specific  boundary  geometries  be  used  in  model  studies  for  seepage 

anal ysi  s. 
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CHAPTER 


CONCLUSSION 


In  this  day  and  age,  more  and  more  sites  are    chosen  for 
buildings  and  engineered  structures  that  are  less  desirable  to 
work  with  from  a  standpoint  of  site  conditions,  poor  soils  and 
water  problems.   A  less  desirable  site  may  easily  be  where  the 
groundwater  table  is  very  near  to  the  surface.   Therfore  the 
building's  design  must  incorporate  some  thought  into  how  the 
groundwater  will  be  handled  and  what  effects  it  may  have  on  the 
structure. 

Engineers  faced  with  the  analysis  of  seepage  for  buildings  to 
be  built  below  the  groundwater  table  must  consider  such 
conditions  of  uplift  pressures,  critical  velocities,  seepage 
control  of  the  water  itself  or  particle  migration.   Hopefully 
the  design  will  enhance  stability  of  the  structure  and  prolong 
the  usefulness  and  life  of  the  structure. 

In  performing  the  literature  review  on  this  topic,  it  was  found 

that  many  authors  and  researchers  a.re    more  and  more  relying  on 

and  developing  computor  models  to  perform  the  analysis  of 

complex  seepage  problems.   With  the  use  of  general  seepage 

theory  and  mathematical  representation  of  the  seepage  theories, 
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the  computer  programs  can  quickly  and  efficiently  aid  the 
endgineer  in  studying  the  flow  conditions.   As  in  this  study, 
the  analysis  of  one  particular  seepage  problem  using  a  -finite 
element  program  made  it  possible  to  study  many  key  variables 
and  develope  many  relationships. 

For  this  study,  underdrain  systems  -for  buildings  built  below 
the  groundwater  table  were  researched.   By  using  a  -finite 
element  program,  and  by  developing  a  mesh  -for  this  particular 
•flow  matrix,  many  variables  o-f  the  analysis  were  studied.   It 
was  -found  that  the  model  could  easily  be  used  to  efficiently 
and  accurately  study  these  variables.   Of  major  concern,  the 
dicharge  of  the  system  and  location  of  the  free  surface  was 
studied  in  detail  versus  the  various  varible  analysis. 

This  research  has  attempted  to  present  general  seepage  theory 
and  useful  criteria  necessary  to  analyzing  this  particular 
seepage  problem.   It  has  also  attempted  to  use  state  of  the  art 
techniques,  i.e.  finite  element  models,  to  analyze  the 
variables  of  the  underdrain  system  design.   In  conclusion,  the 
author  must  urge  that  in  using  such  models  that  the  basic 
assumptions  of  the  model  be  thoroughly  understood  due  to  the 
somewhat  simplifications  of  the  some  model  operations,  and  that 
the  engineer  obtain  as  much  information  and  data  for  the  soil 
parameters  and  fluid  boundary  conditions  which  may  exist. 
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APPENDIX  A-l 

USER'S  MANUAL  FOR 
ARAL  SEEPAGE  PROGRAM 

I.   INTRODUCTION: 

This  seepage  program  is  a  -finite  element  model  for  solving 
uncon-fined  or  con-fined  two  dimensional  or  axi  symetr  i  c  seepage 
problems.   For  uncon-fined  -flow  analysis  the  model  iterates  by 
moving  the  free  surface  nodes  to  the  phreatic  level.  i 

In  general,  the  program  allows  for  the  porous  soil  mass  to  be 
divided  into  discrete  elements  and  nodes  using  either  a  vertical 
plane  model  (two-dimensional  ,>:-y)  or  axisymmetric  grid  structure. 
Then  applying  known  boundary  conditions,   the  program 
simultaneously  solves  for  each  element,  the  continuity  equation 
which  satisfies  Darcy's  Law. 

Figure  1  shows  a  very  simple  case  of  a  two  dimensional  unconfined 
steady  state  condition  which  can    serve  as  an  example.   In  this 
example  the  free  surface  location  is  unknown  and  the  user  must 
begin  the  process  by  i ni tial 1 izi ng  the  input  by  approximating  the 
location.   Then  the  program  will  iterate  to  the  free  surface  level 
until  the  free  surface  converges  to  satisfy  the  equations  of 
flow. 
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Figure  1   Example  Grid  -for  Uncon-fined  Flow  Continu 
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II.  PROGRAM  EXECUTION  AND  ANALYSIS: 

The  user  has  the  option  of  performing  a  two  dimensional  problem  or 
an  axisymmetric  problem.   The  two  dimensional  case  simply  shown  by 
•figure  1  will  serve  as  an  example  to  explain  the  parameters.   An 
axisymetric  application  might  be  a  well  source  problem  looking  at 
stratified  soils  having  varying  permeabilities  at  different 
elevations  through  the  porous  media. 

III.  BOUNDARY  CONDITIONS  OF  FLOW: 

The  potential,  0,  of  flow  must  be  defined  as  some  directional 
function 


JZJ  =   /<x,y,z,t) 


A  Dirichlet  boundary  condition  is  where  the  piezometric  head,  h  i 
constant,  or 

h  =  X(x,y,z,t) 


/<*, 


is  solved  for  at  a  particular  point.   This  boundary  condition 
occurs  whenever  the  flow  domain  is  adjacent  to  an  open  body  of 
fluid.   Since  the  piezometric  head  is  constant  at  the  open  body  of 

water,  the  piezometric  head  of  the  porous  media  at  the  interface 

A- 3 


will  be  constant.   This  boundary  condition  denotes  a  non-  zero 
flux  or  -flow  is  to  be  considered  possible. 

Another  boundary  condition  -for  flow  potential  is  the  Neuman 
boundaries.   This  is  the  case  where  -flow  cannot  move  through  the 
porous  media  at  a  specific  point  defined,  i.e.  an  impervious  lay 
of  stratum. 

IV.   FINITE  ELEMENT  GRID  NOTATION: 

The  porous  soil  mass  to  be  analysed  must  be  divided  into  a     -: 
grid  of  elements  and  nodes.   The  program  uses  any  shape  second-order 
isoparametric  element,  shown  by  figure  2  as  a  quadral ateral 
element  containing  four  corner  nodes  and  four  midpoint  nodes. 


CORNER  NODES 


MID-NODES 


Figure  2   Element  and  Node  Notation 
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The  nodes  must  be  numbered  -from  some  beginning  corner,  preferably 
the  number  one  node  should  be  at  the  bottom  le-ft  corner  o-f  the 
mesh.   Then  node  numbering  increases  to  the  end  o-f  each  column, 
beginning  at  the  bottom  o-f  each  column  throughout  the  mesh.   One 
can  see  the  numbering  scheme  by  the  examples  presented.   Elements 
are  numbered  the  same  as  the  nodes. 

The  program  automatically  generates  nodes  at  equal  spacing  between 
corner  nodes  which  are    defined  by  the  user.   This  feature      -: 
requires  that  careful  attention  be  made  to  the  node  numbering  scheme 
once  corner  nodes  have  been  established.   This  feature  is  good  for 
creating  finer  mesh  grids  in  particular  areas  of  interest. 

Base  elements  are  used  to  create  stationary  reference  elements 
typically  at  the  bottom  of  the  mesh  or  away  from  the  free  surface 
elements.   These  elements  are  also  used  to  define  changes  in 
hydraulic  conductivity  in  the  x  and  y  di  rect  ions  (Km  and  Ky)  . 

The  free  surface  or  phreatic  surface  of  a  media  is  defined  by  the 
user  as  a  first  approximation  of  it's  actual  location.   In  the  case 
of  a  cutoff  wall,  shown  in  Figure  3,  there  is  no  free  surface  and 
the  seepage  model  solves  for  the  flow  characteristics  of  each 
element  and  node.   However,  in  the  case  of  a  dam  as  shown  in 
Figure  3,  the  free  surface  is  unknown  and  the  program  will  solve 

for  the  phreatic  surface  where  the  fluid  pressure  is  at 
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atmospheric  pressure.   In  con-fined  aqui-fers  the  water  pressures 
would  possible  be  higher  than  atmospheric  pressure. 

Seepage  -faces  are    de-fined  by  the  user  to  denote  a  boundary  or 
where  seepage  is. to  be  allowed.   Here  -flow  is  quhatified  and  the 
model  solves  for  the  velocity  and  discharge  at  this  face. 


CUTOFF  WALL 


DAM    CROSS -SECTION 


Figure    3      Example    of    confined    flow    and    unconfined    flow. 
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Permeability  o-f  the  porous  media  often  depends  on  the  direction 
of  -flow.   For  the  condition  of  soil  matrix  being  anisotropic  the 
velocities  and  gradients  Are    generally  not  parallel.   For  this 
condition  the  velocities  can  be  written  in  the  Darcy  Law  -form: 


Vx 


Vy 


-K>;>:  _ 


-Kyy 


Vz   =   -Kzz 


Normally  the  hydraulic  conductivity  of  teh  media  is  similar  in  the 
y  direction.  Therfore  treatment  of  this  normal  condition  requires 
use  of  only  two  direcetions.  For  this  simplification,  the  program 
uses  the  notation  of  the  two  dimensional  case.  Figure  4  shows  the 
notation  that  should  be  used  for  describing  the  hydraulic 
conductivity  conditions. 


Figure  4    Hydraulic  Conductivity  Notation. 
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V.    INPUT  VARIABLES  AND  DATA  INPUT 

Input  can  be  given  in  un -for mated  notation  where  each  data  input  is 
seperated  by  a  comma  or  space.   Careful  attention,  however,  must 
be  made  to  use  the  apppropriate  number  input  data  per  line  as 
described  by  this  guide. 

Title  and  Type  Problem  Card. 

TITLE    The  title  can  be  up  to  72  characters. 

DIMEN    Flow  problem  e  problem  is  a  two  dimensional  analysis 
input   2  and  -for  an  axisymetric  problem  input  1. 

Plane  Flow  Problem      2 

Axisymetric  Problem     1 

FERR     Having  performed  a  sensitivity  of  the  model  or  decided 
on  an  error  critiria,  input  the  error  required. 

Control  cards 

NNPC     Total  number  of  corner  nodal  points. 

These  nodes  are    used  to  define  corners  of  elements 
which  are    defined  by  input  coordinates.   Nodes  that 
are    not  listed  will  be  generated  between  two  corner 
nodes  for  each   vertical  column  or  row  deliniation. 

NELEMC   Total  number  of  base  elements. 

Base  elements  are    typically  at  the  bottom  of  a  series 
of  elements  that  change  the  hydraulic  conductivity  of 
the  elements  above  these  base  elements.   These  elements 
are    established  as  reference  elements  for  the  program 
to  describe  permeability  of  each  element. 

NPBOC    Total  number  of  Neuman  Boundary  Faces. 

A  Neuman  face  is  defined  as  a  line  of  nodes  that  may 
create  a  barrior  or  impervious  wall  for  the  flow.   If 
two  corner  nodes  define  a  line  of  several  nodes  that 
makes  an  impervious  barrier,  this  is  one  Neuman 
Boundary  face. 

NPFS     Total  number  of  corner  nodes  on  the  free  surface. 

The  user  must  specify  particular  corner  nodes  which 
reflect  the  first  approximation  of  the  unknown  free 
surface.   The  program  uses  these  nodes  and  the  midpoint 
nodes  to  converge  to  the  free  surface. 
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ITGIv"    Number  of  iterations  desired. 

This  is  a  safety  -feature  -for  possible  input  errors  and 
computor  time  limitation. 

NPBDC    Total  number  of  Dirichlet  Boundary  -faces. 

Flow  is  allowed  at  these  faces.   Do  not  count  the  nid 
point  nodes  as  part  o-f  the  -face. 

NSPFAC   Total  number  of  seepage  -faces  with  Dirichlet 

boundaries.   For  each  face  defined  by  two  corner  nc  es 
for  the  NPBDC  input  variable  a  NSPFAC  value  should  be 
given.   As  i n 

NPBDC  =  2     NSPFAC  =  1 

NPBDC  =  3     NSPFAC  =  2 

ITIME   Time  variable.   Input  0  =  time  independent. 

This  is  used  for  transient  computation  and  is  not 
used  in  steady  state  computations. 


Corner  Node  Location  Cards 

Corner  node  number  (n) 

x  (n)   x  coordinate  of  node(n) 

y(n)   y  coordinate  of  node(n) 

Note:   There  should  be  NNPC  number  of  these  cards.   Use 
appropriate  dimensions,  i.e.  feet,  meters  which  will  be  carried 
throughout  the  remaining  input  variables. 

NPMIS    Enter  a  "1"  or    "0"    where 

1  -  additional  nodes  are  to  be  generated  after  this 

node. 
0  -  no  nodes  will  be  generated  after  this  node. 
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Base  Reference  Elements  Identification  Cards 


NOD 


For  each  base  element (N),  give  a  counter-clockwise 
numbering  tor  each  node  in  the  element. 


3 

13 

20 

'2 

i 

'19 

1 

> 

i?    , 

18 

NMIS     Total  number  o-f  elements  above  the  base  element  before 
the  next  base  element. 

XK(N)    Hydraulic  conductivity  of  the  elements  defined  by  the 
base  element  in  the  x  direction. 

YK(N)    Hydraulic  conductivity  in  the  y  direction. 

ALP(N)   Angle  of  possible  stratification  or  plane  of  soil 
matrix  with  the  horizontal. 


*-x 


Note:   There  should  be  NELEMC  number  of  these  cards  provided. 


Neuman  Boundary  Condition  Cards 

NSTART   The  beginning  or  lowest  corner  node  number  of  each 
Neuman  Boundary  defined. 

BVAL     The  Neuman  Boundary  value  that  defines  the  known 
potential  height  or  water  level  measured  from  the 
datum  reference  line  of  the  coordinate  system. 

NBSAME   Total  number  of  boundary  faces  having  same  Neuman 

boundary  condition.   Each  face  is  a  line  made  by  two 
corner  nodes  of  the  same  Neuman  boundary. 
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NINC     The  number  of  nodal  number  increases  or  numerical 

di  f  f  erence  to  the  next  corner  node  -from  node  NSTART. 
Do  not  count   mid-nodes  only  corner  nodes. 

Note:   There  should  be  NPBOC  cards  provided  -for  this  set. 


Dirichlet  Boundary  Condition  Cards 

NSTAR    The  beginning  or  lowest  corner  node  number  which 
defines  each  Dirichlet  boundary. 

DBVAL   The  Dirichlet  boundary  value  that  de-fines  the  potential 
or  water  level  height  measured  -from  the  data  reference 
elevation  o-f  the  coordinate  system. 

NDSAME   Total  number  of  boundary  faces  of  this  Dirichlet 
boundary  condition.   Each  face  is  counted  between 
corner  nodes  defining  the  boundary  or  number  of 
elements. 

NDINC    The  number  of  nodal  number  increases  or  numerical 

difference  to  the  next  corner  node  from  node  NSTAR. 

NMID     The  number  of  nodal  number  increases  from  the  NSTAR  to 
the  next  mid-node. 

Note:   There  should  be  NPBDC  cards  provided  for  this  set. 

Free  Surface  Nodal  Cards  (Maximum  16  pairs  per  card) 

NFFSA    Free  surface  corner  node  number. 

NPBOT    Node  number  of  the  bottom  of  the  movable  section  of 
the  grid  that  defines  the  limits  of  moving  nodes  for 
NFFSA. 

Note:   There  should  be  NPFS  pairs  of  node  numbers  given  for  a 
maximum  of  16  pairs  per  card.   Use  as  many  cards  as  required.   The 
bottom  node  (NPBOT)  defines  the  bottom  of  a  seepage  face 
corresponding  to  the  grid  geometry  and  corner  node  defining  the 
free  surface.   The  program  uses  this  band  of  nodes  and  elements  to 
converge  on  the  free  surface. 
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Mid-point  Free  Surface  Nodal  Cards. 

NPFSM    The  mid-point  nodal  numbers  between  corner  nodes 
de-fining  the  -free  surface  (NPFS  -  1). 


Seepage  Face  Cards 

NSPC     Beginning  node  number  of  the  defined  seepage  face. 

NSSAME   Number  of  seepage  faces  defined  by  elements  for  this 
seepage  face  NSPC. 

NSINC    Node  numbering  increase  between  nodes  defined  by 
NSSAME. 

Note:   There  should  be  NSPFAC  cards  provided  for  this  set. 

Statement  and  End  Card  - 

End  of  Problem  Statement  should  be  less  than  72  characters  and 
end  with  a  period. 

NDIMEN   Input  a  "3"  to  end  the  problem. 

ERR      Specify  any  error  as  for  FERR. 
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UNCONFINED  FLOW  EXAMPLE  PROBLEM 

J H 


0Ase  x&eex&jce 


-SA?/ 


69,  79,  86  and 


Title  and  Type  Problem  Card 

TITLE  -  SILTY  SAND  DAM 

DIMEN  -  2  (two  dimensional) 

FERR   -  0.01  (error  criterion  is  0.01) 

Control  Cards 

NMPC   -  12  (corner  nodes  Are    18,  52  35,  45,  52,  6Z 

96) 
NELEMC-  5  (  -five  base  elements  -  shaded) 
NPBOC  -  0  (  no  Neuman  Boundaries) 
NPFS   -6(6  nodal  points  on  the  -free  surface  -  11,  28,  45,  62, 

79  and  96) 
ITGIV  -  50  (  maximum  number  o-f  iterations) 
NPBDC  -  2  (There  Are    two  faces  with  Dirichlet  boundary  conditions, 

the  left  side  and  the  right  side  of  the  mesh  which 

forms  the  dam) 
NSPFAC-  1  (There  is  one  seepage  face  -  on  the  right  side  where 

flow  is  to  be  considered.) 
ITIME  -  0  (Time  independent  for  steady  state  conditions) 
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Corner  Node  Location  Cards 

1,0, 0 ,  1 

Where,  the  -first  node  is  number  1  having  an  «  coordinate  of  0  and 
a  y  coordinate  o-f  0.   Generation  of  additional  nodes  will  -follow. 

11,0, 110, 0 

Where,  node  11  i s  at  the  end  of  the  first  column  at 
coordinates  0  and  110  and  no  generation  is  required  for  other 
nodes  in  this  column. 


18, 200, 0, 1 


28, 

200 

,99,0 

■-•J, 

400 

,  0,  1 

45, 

400 

,  38 ,  0 

52, 

600 

.0,1 

62, 

600 

,77,0 

69, 

800 

,  0,  1 

79, 

800 

,  66 ,  0 

86, 

1 000 ,0,1 

96, 

100< 

J,  _>._>,  U 

Base  Reference  Elements  Identification  Cards 

1, 12, 18, 19,20, 13,3,2,4,0.002.0.002,0 

Counter  clockwise  numbering  of  nodes  for  each  element  and  the 
permeability,  0.002  ft/day  for  this  base  element  and  4  elements 
above  the  base  element. 


1 8 , 29 , 35 , 36 , 37 , 30 , 20 , 1 9 , 4 , 0 . 002 , 0 . 002 , 0 
35 , 46 , 52 , 53 , 54 , 47 , 37 , 36 , 4 , 0 . 002 , 0 . 002 , 0 
52,63,69, 70, 71,64,54,53,4, 0.002,0.002,0 
69,  80,  86,  37,  88,  81 ,  71  ,  70,  4,  0.  002,  0.  002,  0 

Neuman  Boundary  Condition  Cards 

None  specified,  so  no  cards  will  be  provided 
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Dirichlet  Boundary  Condition  Cards 

1, 110,5,2, 1 

86 ,  jj, j, 2,  1 

The  -first  node  of  the  first  Dirichlet  Boundary  is  node  1  and  the 
DBVAL  value  is  110  feet.   There  are    5  elements  on  this  -face    that 
makes  this  Dirichlet  boundary  condition,  NDSAME  =  5.   NDINC  equals 
2  for  the  nodal  number  increase  -for  the  consecutive  elements,  as 
in  3  -  1,  5-3,  7-5,  etc.   NMID  equals  1  for  the  node  numbering 
increase  between  successive  nodes. 

Free  Surface  Nodal  Cards 

11, 1,28, 18,45,35,62,52,79,69,96,86 

Here  the  total  mesh  is  allowed  to  move,  where  the  user  may  -fix  the 
bottom  nodes  at  any  level  in  the  mesh. 

Mid-point  Free  Surface  Nodal  Cards 

17,34,51,68,85 

Seepage  Face  Cards 

86,2,2 

The  seepage  is  prescribed  to  start  at  node  86,  NSFC  =  86,  and 
there  are    2  elements  for  this  seepage  face.   NSINC  is  equal  to  2 
as  there  is  a  two  nodal  increase  between  element  corner  nodes,  as 
is  88  -  86  and  90  -  88. 

Statement  and  End  Card 

END  OF  SILTY  SAND  DAM. 
3 , . 00 1 
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SILTY  SAND  DAM  . 

2,0.01 

12,5,0,6,50,2,1,0 

1,0,0,1 

11,0,110,0 

18,200,0,1 

28,200,99,0 

75,400,0,1 

*5,400,88,0 

62,600,0,1 

62,600,77,0 

fj9,800,0,l 

79,800,66,0 

B6, 1000,0,1 

96,1000,55,0 

1,12,18,19,20,13,3,2,4,0.002,0.002,0 

18,29,35,36,37,30,20,19,4,0.002,0.002,0 

35,46,52,53,54,47,37,36,4,0.002,0.002,0 

52,63,69,70,71,64,54,53,4,0.002,0.002,0 

59,80,86,87,88,81,71,70,4,0.002,0.002,0 

1,110,5,2,1 

36,55,5,2,1 

11,1,28,18,45,35,62,52,79,69,96,86 

17,34,51,68,85 

36,2,2 

2ND  OF  SILTY  SAND  DAM. 

3,0.01 
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iIMCN   ■ 


S1LTY  SAND  OAM 

2 


NNPC  •   12     NELCMC  ■    S     NPBOC   s    o     NPFS  ■    6     ITCIV  '  SO 
GENERATED  NODAL  POINT  OATA  AND  ( X ) - ( V )  COORDINATES 


NSPFAC  • 


NODE 

1 

x  ■ 

O 

OOOO 

Y  i 

0 

OOOO 

NODE 

3 

X  f 

0 

0000 

¥  s 

22 

OOOO 

NOOE 

S 

X  i 

0 

0000 

Vi 

44 

OOOO 

NODE 

7 

X  • 

0 

OOOO 

V  i 

SS 

OOOO 

NOOE 

1 

x  * 

0 

OOOO 

V  I 

ss 

OOOO 

NODE 

1  t 

X  i 

0 

OOOO 

y  i 

1  10 

OOOO 

NODE 

It 

X  ■ 

200 

OOOO 

Yi 

0 

OOOO 

NODE 

20 

X  ■ 

200 

OOOO 

Yi 

It 

tooo 

NODE 

22 

x  ■ 

200 

OOOO 

Y  t 

39 

tooo 

NODE 

24 

X  = 

200 

OOOO 

Yi 

59 

4000 

NOOE 

26 

X  i 

200 

OOOO 

Y  i 

79 

2000 

NODE 

28 

X  t 

200 

OOOO 

Y  • 

St 

OOOO 

NOOE 

35 

X  t 

400 

OOOO 

Yi 

0 

OOOO 

NODE 

37 

X  i 

400 

OOOO 

Y  i 

17 

tooo 

NOOE 

39 

Xi 

400 

OOOO 

Yi 

35 

2O00 

NODE 

4  1 

X  t 

400 

OOOO 

Yi 

52 

tooo 

NOOE 

43 

X  i 

400 

OOOO 

Yi 

70 

4000 

NODE 

4S 

X  ■ 

400 

OOOO 

Yi 

St 

OOOO 

NODE 

S2 

x< 

SOO 

OOOO 

Y  i 

0 

OOOO 

NODE 

54 

Ki 

GOO 

OOOO 

Yi 

IS 

4000 

NOOE 

56 

X: 

600 

OOOO 

Y  I 

30 

sooo 

NOOE 

se 

X  I 

SOO 

oooo 

Y  I 

4G 

20O0 

NODE 

•  0 

X  = 

SOO 

0000 

Y  I 

S  1 

sooo 

NODE 

S2 

X  : 

SOO 

OOOO 

Y  I 

77 

OOOO 

NOOE 

•  9 

X  = 

SOO 

OOOO 

Y  i 

0 

OOOO 

NOOE 

7  1 

X  ■ 

800 

OOOO 

Y  i 

1  3 

2000 

NOOE 

73 

X  I 

800 

OOOO 

Y  i 

26 

4000 

NODE 

75 

X  • 

800 

OOOO 

Y  i 

39 

sooo 

NODE 

77 

X  • 

800 

OOOO 

Y  i 

52 

sooo 

NOOE 

7S 

X  > 

SOO 

OOOO 

Yi 

SS 

OOOO 

NODE 

t* 

x> 

lOOO 

OOOO 

Y  • 

0 

OOOO 

NOOE 

•  S 

Xi 

lOOO 

OOOO 

Y  I 

1  1 

OOOO 

NODE 

SO 

Xi 

1000 

OOOO 

Y  t 

22 

OOOO 

NOOE 

S2 

X  I 

1000 

OOOO 

Y  I 

33 

OOOO 

NODE 

S4 

Xi 

1000 

OOOO 

Yi 

44 

OOOO 

NODE 

St 

Xi 

lOOO 

OOOO 

Y  I 

SS 

OOOO 

NUMBER  OP  NODAL  POINTS...    96 

NUMBER  OF  ELEMENTS 26 

ELE.  NO         NOOAL  POINTS 


1 

12 

It 

1  t 

20 

13 

3 

2 

2 

1  3 

20 

21 

22 

14 

S 

3 

14 

22 

23 

24 

15 

7 

4 

15 

24 

25 

26 

1  S 

t 

5 

1  6 

26 

27 

28 

17 

1  1 

10 

6 

1  S 

29 

35 

31 

37 

30 

20 

1  t 

7 

20 

30 

37 

31 

39 

31 

22 

21 

S 

22 

31 

31 

40 

4  1 

32 

24 

23 

9 

24 

32 

4  1 

42 

43 

33 

26 

25 

10 

26 

33 

43 

44 

49 

34 

21 

27 

1  1 

35 

4f 

52 

53 

54 

47 

37 

38 

1  2 

37 

47 

54 

55 

55 

48 

39 

31 

13 

39 

4t 

St 

57 

St 

41 

41 

40 

14 

41 

4t 

51 

St 

to 

to 

43 

42 

15 

43 

SO 

to 

91 

82 

51 

45 

44 

1  1 

52 

13 

tt 

70 

71 

54 

54 

S3 

17 

54 

t4 

71 

72 

73 

II 

56 

55 

It 

55 

tl 

73 

74 

75 

tt 

51 

57 

It 

SS 

tl 

79 

79 

77 

97 

to 

St 

20 

SO 

17 

77 

71 

79 

tt 

62 

6  1 

21 

tt 

to 

tt 

17 

tt 

t  1 

71 

70 

22 

71 

t  1 

tt 

tt 

to 

12 

73 

72 

23 

73 

82 

to 

9  1 

92 

13 

75 

74 

24 

75 

83 

92 

93 

94 

•  4 

77 

7S 

25 

77 

84 

tt 

ts 

99 

•  5 

75 

71 

OIRICHLET  50UNDARY  CONDITION  DATA 
OES  RESPECTIVE  OIRICHLET  BOUNDARY  CONDITIONS 


2 

1 10 .0000 

1 10  OOOO 

4 

1 10. OOOO 

1 10. OOOO 

S 

110. OOOO 

1 10  OOOO 

1 

110. OOOO 

110. OOOO 

10 

1 10 .0000 

1 10.0000 

10 

1 10.0000 

1 10  OOOO 

97 

55 . OOOO 

55  OOOO 

99 

55 . OOOO 

55  OOOO 

•  91 

55 . OOOO 

55  OOOO 

93 

55 .0000 

55 . OOOO 

95 

SB .OOOO 

55 . OOOO 

99 

55  OOOO 

55  OOOO 

. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 2OOO0E -02 
. 20000E-02 
. 20000E-02 
. 2000OE-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20000E-02 
. 20O00E-02 
. 20O00E-02 
. 20000E-02 


. 20000E 

■02 

0 

. 20OOOE 

■02 

O 

. 20000E 

■02 

O 

. 20000E 

■02 

O 

. 20000E 

■02 

O 

. 20OOOE 

■02 

0 

. 20000E 

■02 

o 

. 20000E 

■02 

o 

. 200OOE 

■02 

o 

. 20000E 

■02 

o 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 2O0OOE 

•02 

0 

. 2O0OOE 

■02 

0 

. 2000OE 

■02 

0 

. 20000E 

■02 

o 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

. 20000E 

■02 

0 

>»L  POINTS  ON  FREE  SURFACE 
II     17     26     34     45     51 


SEEPAGE  FACE  NODAL  POINT  OATA 
INSEOUTIVE  NOOE  NUMBERS  ON  THE  SEEPAGE  FACE 


91 

10 
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13 
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0 
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OOOO 

CONVERGENCE  CHECK  ERROR 


1 7  I OE -02 


X  COORDINATE 


Y  COORDINATE 
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oooo 

1  00 

oooo 

1  00 
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oooo 

400 
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SOO 
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BOO 
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SOO 
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soo 
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SOO 
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SOO 

oooo 
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oooo 

BOO 
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o 
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22 
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8E 
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21 

1  7  1  S 

42 
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63 
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ESEO 

05 

8S7S 

0 

OOOO 

1  0 
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20 

3430 

30 
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50 
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B4 

BE  1  4 

73 
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E  -03 
E  -03 
E  -03 
E-03 
E  -03 
E  -03 


it          600 

OOOO 

1  6 

3SS4 

II           BOO 

OOOO 

24 

5330 

>•           SOO 

oooo 

32 

7  107 

17           BOO 

oooo 

to 

6684 

11           600 

oooo 

49 

06S  1 

11          6  00 

oooo 

57 

2436 

10            BOO 

oooo 

66 

42  14 

1  1           BOO 

oooo 

73 

S99  1 

12           600 

oooo 

6  1 

77BB 

13           700 

oooo 

0 

OOOO 

14           7O0 

oooo 

1  s 

13  14 

IS           700 

oooo 

30 

2628 

16           70O 

oooo 

45 

3942 

17           700 

oooo 

•  0 

52S6 

iS           700 

oooo 

75 

6S70 

jit           BOO 

oooo 

0 

OOOO 

'0          BOO 

oooo 

6 

9537 

M        boo 

oooo 

1  3 

9074 

'2           800 

oooo 

20 

86  11 

3          BOO 

oooo 

27 

8  1  49 

4           BOO 

oooo 

34 

7686 

5          BOO 

oooo 

4  1 

7223 

'6          BOO 

oooo 

48 

67SO 

7           800 

oooo 

5  5 

S297 

6          800 

oooo 

62 

5834 

9           SOO 

oooo 

69 

537  1 

0            900 

oooo 

0 

OOOO 

900 

oooo 

1  2 

4537 

2           900 

oooo 

24 

9074 

3          900 

oooo 

37 

36  11 

4           900 

oooo 

49 

8  149 

S          900 

oooo 

62 

2666 

1         1  000 

oooo 

0 

OOOO 

7         1  OOO 

oooo 

5 

SOOO 

1         1  OOO 

oooo 

1  1 

OOOO 

9         1  000 

oooo 

1  6 

SOOO 

0         1  OOO 

oooo 

22 

OOOO 

1          1  000 

oooo 

27 

5000 

2         1  000 

oooo 

33 

oooo 

3         1  OOO 

oooo 

38 

sooo 

4         1  000 

oooo 

44 

oooo 

5         lOOO 

oooo 

49 

sooo 

6         1  OOO 

oooo 

55 

oooo 

.6166 
.  6234 
632S 
.644  1 
.6584 
.  67S  1 
6944 
.  7  1  80 
.  7401 
75  8434 
75  .  6489 
75 . S6SS 
75 . 6932 
75 . 9322 
75 . 9827 
69  .  6922 
69 . S93S 
69 . 5980 
B9  6049 
69.6151 
69  627S 
69  S43 1 
69 . 6606 
69.68  10 
69 . 7030 
69 . 7274 
62 . 88S6 
B2 . 692S 
62.7129 
62 . 7475 
62 . 7972 
62 . 6634 
SS . OOOO 
55 . OOOO 
55 . OOOO 
SS . OOOO 
55 . OOOO 
SS . OOOO 

ss . oooo 

55 . OOOO 
55 . OOOO 
5  5 . OOOO 
5  5 . OOOO 


1  10SE 

•03 

1  106E 

■03 

1  10SE 

•03 

1  lOSE 

■03 

1  104E 

■03 

1  103E 

•03 

1  102E 

■  03 

1  101E 

•03 

1  1  OOE 

•03 

1202E 

•03 

1  202E 

■03 

1  20  1  E 

■03 

1  200E 

03 

1  1  99E 

■03 

1  1  97E 

•03 

1  301  E 

•03 

1  301E 

•03 

1  300E 

•03 

1  299E 

■03 

1  298E 

03 

1  29SE 

■03 

1  294E 

■03 

1  292E 

■03 

1  288E 

■03 

1  286E 

•03 

1  282E 

03 

1  4S9E 

03 

1  4EOE 

03 

1  4E0E 

03 

1  4E2E 

•03 

1  4E4E 

•03 

1  4E6E 

03 

1  E  1  SE 

03 

1  B  1  SE 

•  03 

1  6  1  7E 

03 

1  620E 

■  03 

1  E23E 

■03 

1  E29E 

■03 

1  E34E 

■03 

1  E42E 

03 

1  649E 

03 

1  EE  1  E 

03 

1  ESSE 

03 

I  2 

I  9 

.  25 

3  1 
38 

4  4 
50 
SS 
62 
32 
1  4 
29 
43 
57 
72 
1  6 
82 
I  8 
24 
36 
40 
54 
54 
70 
66 
66 

.  1  1 

,  1  7 

34 

52 

7  1 

9  4 

.  17 

,  10 

50 

1  0 

1  o 


70E -OS 
IOE -OS 
39E-OS 
7 1 E-OS 
02E -OS 
OSE -OS 
S3E-0S 
S9E -OS 
95E -OS 
1 7E-09 
E 1 E-OS 
1  BE -OS 
58E -05 
82E -05 
33E -OS 
1 SE -07 
77E -06 
49E -OS 
6 1 E-OS 
7 1 E -OS 
24E-OS 
3 1 E -05 
SBE -OS 
83E-05 
70E -OS 
70E -OS 
OBE -OE 
23E -05 
7BE -OS 
90E -OS 
9  BE -05 
37E -OS 
87E -OS 
S8E  -  1  5 
2 1 E-OE 
5BE  -  IS 
34E -OS 

.27E-05 

I  19E -OS 

I2EE -OS 


1  10EE 

-03 

1  1  OBE 

-03 

1  lOBE 

-OS 

1  1  OSE 

-03 

1  10SE 

-03 

1  104E 

-03 

1  103E 

■03 

1  103E 

-03 

1  1  02E 

•03 

1  202E 

-03 

1202E 

■03 

1  201  E 

-03 

1  201E 

-03 

1  200E 

■03 

1  1  99E 

-03 

1301  E 

■03 

1  301E 

-03 

1300E 

■OS 

1299E 

■03 

1  29BE 

■03 

1  297E 

■03 

1  29SE 

■03 

1293E 

-03 

1  290E 

■03 

1287E 

■03 

1  28SE 

■03 

1459E 

■03 

1  460E 

■03 

1  46  IE 

■03 

1463E 

-03 

1  4BEE 

■03 

14E9E 

-03 

1  6  1  SE 

■03 

1  6  1  6E 

■03 

16  176 

■03 

1  620E 

•03 

1  623E 

■03 

1  629E 

■03 

1  634E 

■03 

1  B42E 

■03 

1  650E 

•03 

1  66  1  E 

■03 

1  ESSE 

•03 

SEEPAGE  FACE  COMPUTATION  RESULTS 


AVE .   VEL 


.16  1 60E -03 


1 1000E*02 


1 7776E-02 


TOTAL  DISCHARGE 
88     NE   =     90 


. 1 7776E -02 
VEL  •    . 1 6202E-O3     OS 


1  1  OOOE-0O2 


17622E-02 


TOTAL  DISCHARGE 


35598E -02 


END  OF  SILTV  SAND  DAM. 
OIMEN  >  3  PB«R  • 
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APPENDIX  B-l 

SAMPLE  INPUT 
AND  OUTPUT 


RUN  ADS2          SP  =  25 

FEET 

2,0.0001 

?3, 30, 0,30, 20, 6, 5,0 

1,-340,0,1 

5,-340,34,1 

11,-340,78.0,0 

18,-200,0,1 

22,-200,34,1 

28,-200,78.0,0 

35,-80,0,1 

59,-80,34,1 

!»5,-80,77.0,0 

52,0,0,1 

56,0,34,1 

52,0,74.0,0 

59,50,0,1 

73,50,34,1 

79,50,71.0,0 

J6, 58. 75, 0,1 

)0, 58. 75, 34,1 

>6, 59. 0,58. 0,0 

103,62.0,0,1 

107,62.0,34,1 

113,61.0,58.0,0 

120,66.0,0,1 

124,66.0,34,1 

130,65.0,60.0,0 

137,70.0,0,1 

141,70.0,34,1 

147,69.0,60.0,0 

154,74.0,0,1 

158,74.0,34,1 

164,75.0,60.0,0 

171,78.0,0,1 

175,78.0,34,1 

181,79.0,60.0,0 

188,82.0,0,1 

192,82.0,34,1 

198,83.0,58.0,0 

'05,86.0,0,1 

109,86.0,34,1 

515,85.0,58.0,0 

122,90.0,0,1 

526,90.0,34,1 

532,89.0,60.0,0 

539,94.0,0,1 

543,94.0,34,1 

549,93.0,60.0,0 

556,98.0,0,1 

560,98.0,34,1 

566,99.0,60.0,0 

573,102.0,0,1 

577,102.0,34,1 

583,103.0,60.0,0 

'.90,106.0,0,1 

594,106.0,34,1 

100,107.0,58.0,0 

^07,110.0,0,1 

HI, 110.0,34,1 

117,109.0,58.0,0 

DEPTH  =   20  FEET     K  =    .2.      DEPTH  =58'. 


B-l 


,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,58.0,0 

,0,1 

,34,1 

,58.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,60.0,0 

,0,1 

,34,1 

,58.0,0 

,0,1 

,34,1 

,58.0,0 

,12,18,19,20,13,3,2,4,0.2,0.2,0.0 
8,29,35,36,37,30,20,19,4,0.2,0.2,0.0 
5,46,52,53,54,47,37,36,4,0.2,0.2,0.0 
2,63,69,70,71,64,54,53,4,0.2,0.2,0.0 
9,80,86,87,88,81,71,70,4,0.2,0.2,0.0 
6,97,103,104,105,98,88,87,4,0.2,0.2,0.0 
03,114,120,121,122,115,105,104,4,0.2,0.2 
20,131,137,138,139,132,122,121,4,0.2,0.2 
37,148,154,155,156,149,139,138,4,0.2,0.2 
54,165,171,172,173,166,156,155,4,0.2,0.2 
71,182,188,189,190,183,173,172,4,0.2,0.2 
88,199,205,206,207,200,190,189,4,0.2,0.2 
05,216,222,223,224,217,207,206,4,0.2,0.2 
22,233,239,240,241,234,224,223,4,0.2,0.2 
39,250,256,257,258,251,241,240,4,0.2,0.2 
56,267,273,274,275,268,258,257,4,0.2,0.2 
73,284,290,291,292,285,275,274,4,0.2,0.2 
90,301,307,308,309,302,292,291,4,0.2,0.2 
07,318,324,325,326,319,309,308,4,0.2,0.2 
24,335,341,342,343,336,326,325,4,0.2,0.2 
41 , 352 , 358 , 359 , 360 , 353 , 343 , 342 , 4 , 0 . 2 , 0 . 2 
58,369,375,376,377,370,360,359,4,0.2,0.2 
75,386,392,393,394,387,377,376,4,0.2,0.2 
92,403,409,410,411,404,394,393,4,0.2,0.2 


124 

114.0, 

128 

114.0, 

134 

113.0, 

141 

118.0, 

145 

118.0, 

151 

117.0, 

158 

122.0, 

162 

122.0, 

168 

123.0, 

175 

126.0, 

179 

126.0, 

(85 

127.0, 

192 

130.0, 

196 

130.0, 

►02 

131.0, 

►09 

134.0, 

►  13 

134.0, 

►  19 

133.0, 

►26 

138.0, 

►30 

138.0, 

►36 

137.0, 

43 

142.0, 

•47 

142.0, 

53, 

141.0, 

60 

146.0, 

64 

146.0, 

70 

147.0, 

77 

150.0, 

81 

,150.0, 

•87 

,151.0, 

94 

,154.0, 

98 

,154.0, 

04 

,155.0, 

11 

,156.0, 

15 

,156.0, 

21 

,156.0, 

,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
,0.0 
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409,420,426 

426,437,443 

443,454,460 

460,471,477 

477,488,494 

494,505,511 

1,78.0,5,2, 

96,58.0,1,1 

198,58.0,1, 

300,58.0,1, 

402,58.0,1, 

504,58.0,1, 

11,5,28,22, 

147,141,164 

283,277,300 

419,413,436 

17,34,51,68 

289,306,323 

96,1,17,198 

SP-25..K  - 

3,0.0001 


,427,428,421,411,410,4,0.2,0.2,0.0 

,444,445,438,428,427,4,0.2,0.2,0.0 

,461,462,455,445,444,4,0.2,0.2,0.0 

,478,479,472,462,461,4,0.2,0.2,0.0 

,495,496,489,479,478,4,0.2,0.2,0.0 

,512,513,506,496,495,4,0.2,0.2,0.0 

1 

7,6 

17,6 

17,6 

17,6 

17,6 

45,39,62,56,79,73,96,90,113,107,130,124 

, 158 , 181 , 175 , 198 ,192,215, 209 , 232 , 226 , 249 , 243 , 266 , 260 

,294,317,311,334,328,351,345,368,362,385,329,402,396 

,430,453,447,470,464,487,481,504,498 

, 85 , 102 , 1 19 , 1 36 , 153 , 1 70 , 187 , 204 , 22 1 , 238 , 255 , 272 

, 340 , 357 , 374 , 39 1 , 408 , 425 , 442 , 459 , 476 , 493 

,1,17,300,1,17,402,1,17,504,1,17 

.2..D  -  20  FT... END  OF  PROBLEM. 


B-3 


1 

RUN  A0S2 

SP   >   25   FEET 

NDIMEN  i 

2 

FERR  > 

DEPTH  s  20  FEET   K  • 
. 000  10 


DEPTH  >  68 


NNPC  •  S3 


NELEMC  <  30 


ITCIV  <  20 


NSPFAC  > 


CENERATED  NODAL  POINT  DATA  AND  IX)-(V)  COORDINATES 


NODE 

1 

Xi 

•340 . 0000 

Yi 

0.0000 

NODE 

3 

»" 

-340 . 0000 

¥• 

17 .0000 

NODE 

I 

X> 

-34O.OO00 

Yi 

34 .OOOO 

NODE 

7 

Xt 

-  340 . OOOO 

Y  ■ 

46 . 6667 

NODE 

9 

Xi 

-340 . 0000 

Yi 

63 . 3333 

NODE 

1  1 

Xi 

-340. OOOO 

Yi 

76 . OOOO 

NODE 

It 

Xi 

-  200 . OOOO 

Yi 

0 . OOOO 

NODE 

20 

Xi 

•200. OOOO 

Yi 

17 . OOOO 

NODE 

22 

Xi 

-200 . OOOO 

Yi 

34 . OOOO 

NODE 

24 

Xi 

-200 . OOOO 

Vi 

44 . 6687 

NODE 

26 

Xi 

-200 . OOOO 

Yi 

63 . 3333 

NODE 

2« 

Xi 

-200 . OOOO 

Y  I 

78  OOOO 

NODE 

35 

Xi 

-SO . OOOO 

Yi 

0.0000 

NODE 

37 

Xi 

-60 . OOOO 

Y  • 

1 7  OOOO 

NODE 

31 

Xt 

-80 . OOOO 

Y  • 

34 . OOOO 

NODE 

41 

Xi 

-60 . OOOO 

Yi 

44 . 3333 

NOOE 

43 

Xi 

-80. OOOO 

Yi 

62 . 8667 

NODE 

45 

Xi 

-80.0000 

Y  • 

77 . OOOO 

NOOE 

82 

Xi 

0 .OOOO 

Y  i 

O . OOOO 

NODE 

54 

Xi 

0 . OOOO 

Yi 

17 . OOOO 

NOOE 

S« 

Xi 

0 . OOOO 

Yi 

34 . OOOO 

NOOE 

5! 

Xi 

0 . OOOO 

Yi 

47 . 3333 

NODE 

60 

Xi 

0 . OOOO 

Yi 

60. 6667 

NDDE 

62 

Xi 

0 .OOOO 

Yi 

74 .OOOO 

NOOE 

69 

Xi 

50 . OOOO 

Y  • 

0 . OOOO 

NODE 

71 

Xi 

SO . OOOO 

Y  • 

1 7 . OOOO 

NODE 

73 

Xi 

SO . OOOO 

Y  • 

34 . OOOO 

NOOE 

76 

Xi 

50 . OOOO 

Yi 

46 . 3333 

NOOE 

77 

Xi 

50 . OOOO 

Yi 

54  6667 

NOOE 

7i 

Xi 

50 . OOOO 

Y  • 

7  1 . OOOO 

NODE 

66 

Xi 

56 . 7500 

Y  i 

0 . OOOO 

NOOE 

66 

Xi 

56 . 7SOO 

Y  i 

1 7 . OOOO 

NOOE 

60 

Xi 

56 . 7600 

Y  • 

34 . OOOO 

NODE 

92 

Xi 

56 . 6333 

Y  ■ 

42 .OOOO 

NODE 

94 

Xi 

58 . 8 167 

Y  I 

SO . OOOO 

NODE 

96 

Xi 

58 . OOOO 

Y  • 

58 . OOOO 

NODE 

103 

Xi 

62 .OOOO 

Yi 

0  OOOO 

NODE 

106 

X  i 

62 .OOOO 

Yi 

17 . OOOO 

NODE 

107 

Xi 

62 . OOOO 

Y  I 

34 .OOOO 

NODE 

109 

X  • 

61 . 8687 

Y  I 

42 .OOOO 

NODE 

1  1  1 

Xi 

61 . 3333 

Y  I 

SO. OOOO 

NODE 

1  13 

Xi 

6 1 . OOOO 

Y  i 

SS .OOOO 

NOOE 

1  20 

Xi 

66 . OOOO 

Y  • 

0 . OOOO 

NODE 

122 

Xi 

66 .OOOO 

Y  • 

17 . OOOO 

NODE 

124 

Xi 

66 .OOOO 

Yi 

34 . OOOO 

NODE 

124 

Xi 

•■ . 6687 

Yi 

42 . 6667 

NOOE 

126 

Xi 

81 .3333 

Yi 

SI . 3333 

NODE 

130 

Xi 

65 .OOOO 

Yi 

60 . OOOO 

NOOE 

137 

Xi 

70 . OOOO 

Yi 

O . OOOO 

NODE 

139 

Xi 

70 . OOOO 

Yi 

17 . OOOO 

NODE 

141 

Xi 

70 .OOOO 

Y  • 

34 . OOOO 

NODE 

143 

Xi 

86 . 6667 

Yi 

42 . 6667 

NODE 

141 

Xi 

69 . 3333 

Yi 

SI . 3333 

NODE 

147 

Xi 

68 .OOOO 

Y  • 

60 . OOOO 

NODE 

164 

Xi 

74 .OOOO 

Yi 

0 .OOOO 

NODE 

166 

Xi 

7  4 .OOOO 

Yi 

1 7 . OOOO 

NODE 

156 

Xi 

74 . OOOO 

Yi 

34 .OOOO 

NOOE 

1  60 

Xi 

74 . 3333 

Y  ■ 

42 . 6867 

NOOE 

162 

XI 

74 . 6667 

Y  • 

SI . 3333 

NOOE 

1  64 

Xi 

75 .OOOO 

Y  • 

60 . OOOO 

NODE 

171 

Xi 

76 . OOOO 

Yi 

0 . OOOO 

NODE 

173 

Xi 

78 . OOOO 

Yi 

17 . OOOO 

NOOE 

175 

Xi 

78 . OOOO 

Yi 

34 . OOOO 

NOOE 

177 

Xi 

76.3333 

Yi 

42 . 4667 

NODE 

179 

Xi 

76. 6667 

Yi 

SI . 3333 

NOOE 

141 

Xi 

76 .OOOO 

Yi 

60 . OOOO 

NOOE 

166 

Xi 

82 . OOOO 

Yi 

0 . OOOO 

NOOE 

190 

Xi 

62 .OOOO 

Yi 

17 . OOOO 

NOOE 

192 

Xi 

82 .0000 

Yi 

34 .OOOO 

NODE 

194 

Xi 

62 . 3333 

Y  I 

42 .OOOO 

NODE 

196 

Xi 

82 . 6667 

Yi 

50 .OOOO 

NODE 

186 

Xi 

83 .OOOO 

Yi 

56  OOOO 

NODE 

205 

Xi 

8  6 . OOOO 

Yi 

O.OOOO 

NODE 

207 

Xi 

8  6 . OOOO 

Yi 

17 . OOOO 

NODE 

209 

Xi 

8  6 . OOOO 

Yi 

34 .OOOO 

NODE 

21  1 

X  « 

86 . 6667 

Yi 

42 . OOOO 

NODE 

213 

Xi 

86 . 3333 

Y  I 

50 . OOOO 

NODE 

215 

Xi 

65  OOOO 

Y  ■ 

S6 . OOOO 

NODE 

222 

Xi 

80 .OOOO 

Yi 

0 . OOOO 

NODE 

224 

Xi 

SO . OOOO 

Yi 

1 7 . OOOO 

NODE 

226 

Xi 

80 . OOOO 

Yi 

34 .OOOO 

NODE 

226 

X  • 

86 . 8887 

Yi 

42 . 6667 

NOOE 

230 

Xi 

68 . 3333 

Yi 

51 . 3333 

NOOE 

232 

Xi 

68 .OOOO 

Yi 

60 . OOOO 

NOOE 

238 

Xi 

84 .OOOO 

Yi 

0 . OOOO 

NOOE 

241 

Xi 

84 . OOOO 

Yi 

17 . OOOO 

NOOE 

243 

Xi 

8  4 .OOOO 

Yi 

34 . OOOO 

NOOE 

246 

Xi 

63 . 6667 

Yi 

42 . 6667 

NODE 

247 

Xi 

63 . 3333 

Yi 

5 1 . 3333 

NODE 

248 

Xi 

S3 .OOOO 

Yi 

60. OOOO 

NOOE 

256 

X  • 

88 .OOOO 

Yi 

O.OOOO 

NODE 

256 

X  . 

86  OOOO 

Yi 

17 . OOOO 

NOOE 

260 

Xi 

8  6  OOOO 

Y  • 

34 . OOOO 

NOOE 

262 

Xi 
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34 
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432 
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137 

6667 
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42 
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434 
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137 

3333 

Y  ■ 
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NOOE 

43S 
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137 

OOOO 
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443 
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142 
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17 
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42 
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141 
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Yi 

10 
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Xi 
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oooo 

Y  i 

0 

OOOO 
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Xi 

1  to 
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Yi 

17 
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Xi 

1  to 
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Yi 

34 
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NOOE 
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Xi 
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3333 

Yi 

42 
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•  85 

Xi 
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Yi 

t  1 

3333 
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Xi 

It  1 

oooo 

Yi 
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Xi 
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Yi 
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Xi 
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oooo 
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Yi 

34 
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114 
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42 
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NODE 
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Xi 

114 
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Yi 

to 
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NODE 

504 

Xi 
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Xi 

156 

oooo 

Yi 

0 

oooo 
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513 

Xi 

155 
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Yi 

17 

.  oooo 

NODE 

515 

Xi 

111 
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Yi 

34 
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NODE 

517 

Xi 

1  tl 
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Xi 
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1 

1 

-  12 

11 

11 

20 

13 

3 

2 

. 20000Et00 

. 200001*00 

O 

2 

3 

13 

20 

21 

22 

14 

S 

4 

. 20000E too 

. 20O00E tOO 

0 

3 

5 

14 

22 

23 

24 

IS 

7 

1 

. 20000E too 

.  20000E  tOO 

0 

4 

7 

15 

24 

25 

21 

IS 

■ 

I 

. 2OO00E too 

. 20OO0E tOO 

0 

1 

1 

11 

21 

27 

21 

17 

1  1 

10 

.  20000E  too 

. 20000E*00 

0 

1 

11 

25 

31 

35 

37 

30 

20 

11 

. 20OO0E too 

. 2OOOOEtO0 

0 

7 

20 

30 

37 

31 

3t 

31 

22 

21 

. 20000E too 

. 20000E+00 

0 

• 

22 

31 

31 

40 

41 

32 

24 

23 

. 20000E*00 

. 20000E+00 

o 

• 

24 

32 

4  1 

42 

43 

33 

21 

25 

. 200001*00 

. 20000Et00 

0 

io 

26 

33 

43 

44 

45 

34 

21 

27 

. 2OO0OE too 

. 20000E tOO 

o 

1  1 

35 

4t 

12 

53 

S4 

47 

37 

35 

. 20000E*00 

. 200O0E too 

0 

12 

37 

47 

S4 

SS 

SS 

41 

31 

36 

. 20000E*00 

. 20000EtOO 

0 

13 

31 

41 

51 

57 

SI 

41 

41 

40 

. 20O0OE too 

. 200006 tOO 

0 

14 

41 

41 

SI 

51 

•  O 

SO 

43 

42 

. 20OO0E too 

. 200OOE tOO 

0 

IS 

43 

to 

to 

■  1 

■  2 

51 

4S 

44 

. 20000E too 

. 20000B too 

0 

It 

12 

13 

II 

70 

71 

64 

S4 

S3 

.  2OO00E t 00 

. 20000E tOO 

0 

17 

54 

•  4 

71 

72 

73 

65 

tl 

55 

. 20000E+00 

. 2OOO0E tOO 

0 

It 

St 

•  t 

73 

74 

75 

51 

tl 

57 

. 20000E*00 

. 2OOO0EtO0 

0 

It 

SI 

11 

71 

71 

77 

•  7 

IO 

SI 

. 20000E*00 

. 20000E tOO 

0 

20 

■  0 

47 

77 

71 

71 

II 

•  2 

11 

. 20000E+00 

. 20000E+00 

0 

21 

It 

to 

•  1 

•  7 

11 

11 

71 

70 

. 20000E*00 

. 20OOOEt00 

o 

22 

71 

•  1 

11 

11 

10 

12 

73 

72 

. 20000E too 

. 20000Et00 

o 

23 

73 

12 

10 

11 

12 

S3 

75 

74 

. 200OOE too 

. 20000E too 

o 

24 

75 

•  3 

12 

13 

14 

•  4 

77 

71 

.  20O0OE  »00 

. 2OO00E too 

0 

25 

77 

•  4 

•  4 

IS 

SI 

It 

71 

71 

. 20000E+00 

.  20000E  too 

0 

21 

•  1 

97 

103 

104 

101 

■  1 

81 

17 

. 200O0E+0O 

. 20000E tOO 

0 

27 

tt 

11 

101 

106 

107 

11 

•  0 

11 

. 20000E*00 

. 20000E tOO 

0 

21 

•  0 

•  1 

107 

101 

101 

100 

•  2 

11 

. 20000E too 

. 20000EtOO 

0 

21 

92 

ioo 

101 

1  10 

1  1  1 

101 

•  4 

13 

. 20000E too 

. 20000E tOO 

o 

30 

•  4 

101 

1  1  1 

1  12 

1  13 

102 

tt 

IS 

. 20000E+00 

. 20000E-00 

0 

31 

103 

1  14 

120 

121 

122 

1  1t 

101 

104 

. 20000E»00 

. 20OOOE too 

0 

32 

10S 

1  It 

122 

123 

124 

1  11 

107 

101 

. 20000E+00 

. 20OO0EtO0 

0 

33 

107 

1  It 

124 

125 

121 

1  17 

101 

101 

. 20OO0E+OO 

.  20000E  tOO 

0 

34 

101 

1  17 

121 

127 

121 

1  11 

1  1  1 

1  10 

. 20OOOE too 

. 20000E tOO 

0 

35 

1  1  1 

1  It 

1  21 

121 

130 

1  11 

1  13 

1  12 

. 20000E too 

. 20000E tOO 

o 

31 

120 

131 

137 

131 

131 

132 

122 

121 

.  20OOOE  too 

. 20000E tOO 

o 

37 

122 

132 

131 

140 

141 

133 

124 

123 

. 20000E too 

. 20000E tOO 

0 

31 

124 

133 

141 

142 

143 

134 

121 

125 

. 20000E+00 

. 20O00E too 

0 

31 

124 

134 

143 

144 

14S 

13S 

121 

127 

. 20000E*00 

. 2O000EtO0 

0 

40 

121 

135 

14S 

141 

147 

136 

130 

126 

. 20000E too 

. 20OO0Et00 

0 

41 

137 

141 

1S4 

155 

158 

141 

131 

138 

. 20000E+00 

.  20O00E  tOO 

0 

42 

131 

141 

156 

1S7 

1SI 

ISO 

141 

140 

. 20000E+00 

. 20000E tOO 

0 

43 

141 

110 

151 

151 

110 

151 

143 

142 

. 20000E+00 

. 20000E tOO 

0 

44 

143 

1  tl 

1  10 

111 

1  12 

152 

US 

144 

. 20000E too 

. 200006 tOO 

0 

4! 

145 

1S2 

1  12 

1  13 

1  14 

153 

147 

146 

. 2OOO0E tOO 

. 20O00Et00 

0 

41 

1S4 

1  IS 

171 

172 

173 

1  66 

151 

1SS 

.  20000E  »00 

. 20000E tOO 

0 

47 

156 

1  11 

173 

174 

17S 

1  67 

158 

157 

. 20000E+00 

. 20000E tOO 

0 

41 

151 

1  17 

17S 

171 

177 

1  61 

1  SO 

151 

. 20000E too 

. 20000E too 

0 

41 

1  10 

1  11 

177 

171 

171 

1  61 

162 

111 

. 20000E too 

. 20000E too 

0 

SO 

1  62 

1  11 

171 

110 

111 

170 

1  64 

1  13 

. 20OOOE+00 

. 20000E too 

o 

SI 

17  1 

112 

111 

111 

ISO 

163 

173 

172 

. 2OO00E+OO 

. 20000E too 

o 

52 

173 

113 

ISO 

111 

1  12 

164 

175 

174 

. 20000E+00 

. 20000E too 

0 

S3 

175 

114 

112 

113 

114 

165 

177 

171 

. 20000E too 

. 20000E too 

0 

S4 

177 

1  IS 

114 

US 

111 

1  11 

171 

171 

.  20000E  too 

. 20OO0E+OO 

o 

55 

171 

1  11 

111 

117 

111 

167 

18  1 

1  10 

. 20000E+00 

. 20000E too 

0 

SS 

1  It 

111 

205 

201 

207 

200 

110 

111 

. 20000E too 

. 2OOO0E too 

o 

S7 

110 

200 

207 

201 

201 

201 

112 

11  1 

. 20000Et0O 

.  2OOO0E  too 

0 

51 

112 

201 

201 

210 

21  1 

202 

114 

113 

.  20000E  tOO 

.  20O0OE  too 

0 

51 

114 

202 

21  1 

212 

21  3 

203 

111 

115 

.  20000E  tOO 

. 20000E too 

o 
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•  4 

•  f 

•  7 

■  • 

■  t 
70 
7  1 
72 
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71 
77 
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71 
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•  1 
12 

•  3 
44 

!  I 

•  t 

•  7 
»• 
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t  1 

■  2 

•  3 
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IS 

•  1 

•  7 
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It 
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lOI 
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104 
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107 
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101 
1  lO 
1  1  1 
1  12 
1  13 
1  14 
1  1* 
1  It 
1  17 
1  II 
1  II 
1  20 
121 
122 
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134 
12* 
121 
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1  21 

:  i 

130 
131 
132 
133 
134 
:  i 
131 
U7 
131 
131 
140 
1  4  1 
142 
143 
1  44 
141 
141 
1  47 
141 
141 

1  to 


207 
20t 
21  1 
213 
222 
224 
221 
221 
230 
231 
241 
243 
241 
247 
2tl 
2tt 
210 
212 
214 
273 
271 
277 
271 
211 
2IO 
212 
214 
211 
211 
307 
301 
31  l 
31  3 
31  t 
324 
321 
324 
330 
332 
341 
343 
34f 
347 
341 
311 
310 
312 
3t4 
III 
371 
377 
371 
31  I 
313 
312 
314 
311 
311 
4O0 
401 
4  1  1 
413 
411 
41  7 
421 
421 


432 
434 
443 
441 
447 
44t 
41  1 
410 
412 
414 
411 
4lt 


217 
2  1  I 
211 
220 
233 
234 
231 
231 
237 
210 
211 
212 
213 
214 
217 
211 
211 
270 
271 
214 
241 
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217 
211 
301 
302 
303 
304 
30t 
311 
311 
320 
321 
322 
331 
331 
337 
334 
331 
312 
313 
3t4 
311 
3tt 
311 
370 
371 
372 
373 
314 
347 
314 
311 
310 
403 
404 
401 
40t 
407 
420 
421 
422 
423 
424 
437 
431 


224 
221 
221 
230 
231 
241 
243 
241 
247 
211 
211 
210 
212 
214 
273 
271 
277 
27t 
211 
210 
212 
214 
211 
2tl 
307 
301 
31  1 
313 

3  l  i 
324 
321 
321 
330 
332 
34  1 
343 
348 
347 
341 
3S4 
:  i  : 
312 
314 
344 
371 
377 
37t 
3tl 
:  i  : 
312 
314 
3tt 
311 
400 
401 

4  1  1 
413 
4  11 
417 
421 
424 
430 
432 
434 
443 
441 


22t 

227 
221 
23  1 
240 
242 
244 
241 
241 
217 
211 
211 
213 
211 
274 
271 
271 
210 
212 
2t  1 
213 
2IS 
217 
211 
301 

3  10 
312 
314 
31  t 
324 
327 
321 
33  1 
333 
342 
344 
344 
344 
3S0 
381 
3f  l 
313 
311 
317 
371 
371 
310 
312 
344 
313 
311 
317 
311 
401 

4  10 
412 
4  1  4 
4  1  I 
I  II 
427 
421 
431 
433 
HI 
444 
44t 


221 
221 
230 
232 
241 
243 
241 
247 
24t 
211 
210 
212 
214 
211 
27S 
277 
271 
211 
213 
212 
2t4 
211 
211 
3O0 
301 
31  1 
31  3 
31  t 
317 
321 
321 
330 
332 
334 
343 
:  I I 
347 
341 
3S1 
340 
312 
314 
311 
311 
377 
371 
311 
313 
3lt 
314 
311 
311 
400 
402 
4  1  1 
413 
411 
4  17 
4  It 
424 
430 
432 
434 
431 
441 
447 


440 
441 
414 
4lf 
411 
417 
411 
471 
472 
473 
474 
471 


441 

411 
410 
412 
414 
411 
411 
177 
471 
411 
413 
411 


410 
412 

411 
413 
411 
417 
411 
471 
410 
412 
414 
411 


41  1 
413 
412 
414 
411 
411 
470 
47t 
411 
413 
•  II 
417 


2  1  I 
211 
220 
221 
234 
231 
234 
237 
231 
2t  l 
212 
:  <  : 

214 
211 
211 
211 

270 
271 
272 
211 
211 
217 
211 
211 
302 
303 
304 
301 
301 
311 
320 
32  1 
322 
323 
331 
337 
331 
331 
340 
313 
314 
3tl 
3tl 
317 
370 
371 
372 
373 
374 
317 
311 
311 
3t0 
31  1 
404 
401 
401 
407 
401 
421 
422 
423 
424 
42t 
431 
43t 


44  1 
442 
411 
411 

417 
411 
411 

472 
473 
474 
471 

471 


201 
21  1 
213 
211 
224 
22t 
221 
230 
232 
241 
243 
241 
247 
241 
211 
210 
212 
214 
2lt 
271 
277 
271 
21  1 
213 
212 
214 
211 
211 
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301 
31  1 
3  1  3 
3  II 
3  17 
321 
321 
330 
332 
334 
343 
341 
347 
341 
31  1 
3IO 

3  1! 
314 
311 
311 
377 
371 
311 
313 
311 
314 
311 
311 
4O0 
402 
•  1  1 

4  13 
411 
417 
4  II 
421 
430 


432 
434 
431 
441 
447 
441 
4(1 
413 
412 
414 
411 
411 
470 


201 
210 
212 

2  14 
223 
22t 
227 
221 
231 
240 
242 
244 
241 
241 
217 
211 
211 
213 
211 
274 
271 
271 
210 
212 
21  1 
213 
211 
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211 
301 

3  10 
3  12 
3  14 
31  t 
321 
327 
32t 
33  1 
333 
342 
344 
341 
344 
310 
311 
3*1 
313 

3  11 
317 
371 
371 
310 
312 
314 
313 
2  11 
317 
311 
401 
410 
412 

4  1  4 
4  1  I 
4  II 
427 
421 
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SEEPAGE  PACE  COMPUTATION  RESULTS 
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TOTAL  DISCHARGE    I         102406*00 
113     NE  >    188     AVE.   VEL  t     381076-01     OS  I 


TOTAL  DISCHARGE    >         940766*00 
186     NE  •    21S     AVE.   VEL  ■    .2S3S2E-01     DS  = 


TOTAL  DISCHARGE    >        .883466+00 
21S     NE  '    300     AVE.   VEL  •     204816-01     DS 


TOTAL  DISCHARGE  144406*01 

300     NE  «    317     AVE.  VEL  >     156496-01     DS 


TOTAL  DISCHARGE    ■         147536*01 
317     NE  t    402     AVE.   VEL  ■     127846-01     DS  ■ 


TOTAL  DISCHARGE    l         175666*01 
402     NE  •    418     AVE.   VEL  >     105636-01     DS  ■ 


TOTAL  DISCHARGE    i         177786+01 
418     NE  ■    S04     AVE.  VEL  •    .88810E-02     OS 


TOTAL  DISCHARGE    >         199516+01 
S04     NE  t    S21     AVE.   VEL  >     90157E-02     DS  ■ 


TOTAL  DISCHARGE  2OO426+01 
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